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SINOPSIS 

This thesis presents a study of proton spin-lattice 

relaxation in pure methane and in mixtures of CH^-N 2 ,CH^-C 02 

and CH^-Inert gas atom (He, He and Ax) using spin-echo technique. 

The measurements were made as a function of density 1 

amagats) and composition in the temperature range 300-600°E, A 

spin-echo spectrometer operating at 30 MHz was fabricated for 

this purpose. The temperature dependences of the spin-lattice 

relaxation time per unit density . ^ T^ , in case of pure CH^ 

and the extrapolated values of (T^/^ ) to 100^ H' 2 » ^2 Inert 

gas atom in case of mixtures along with the data of Armstrong 
18 

and 0?ward in Cl/, and Sil^ have been interpreted to obtain 
information about the intermole cular potentials, especially the 
anisotropic potential in these systems. The analysis is made 
assuming that the average lifetime of a molecule in a J state 
can be approximatod by the correlation time of spin-rotation 
interaction, the dominant mechanism of relaxation in the systems 
under consideration. The average lifetime of a molecule in a 
J state is calculated following the theory of Bloom and Oppenheim 
for a given anisotropic intermolecular potential using ’’weak 
collision approximation''. The molecular quantities such as the 
hyperpolarizability of CH^, the oc to pole moments of CH^, and 
SiP^ have been determined from the experimental data and Ihe 
results are discussed. 
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The first chapter of the thesis introduces the importance 
of gas phase nuclear spin relaxation studies in providing valuabl 
information on the anisotropy of in te molecular interactions. 

The chapter includes an outline of the basic principles of 
nuclear magnetic resonance. Finally the inter and intra mole- 
cular contributions to the relaxation rate are discussed. 

Chapter II deals with the description of the spin-echo 
apparatus fabricated for this work. A brief description of the 
spin-echo technique is given. The essential features of the 
pulsed spectrometer are that the operating frequency is at 30 MHz 
and phase coherent detection is incorporated. Tektronix wave- 
form and pulsewidth generators in combination with a pulse 
sequencer, a pulsewidth- generator of the phantastron type and 
a pulse amplifier produce the necessary sequence of d.c. pulses. 
These d.c. pulses are used to gate a continuous wave crystal 
controlled oscillator at 10 MHz. The gated output is amplified, 
tripled and again power amplified to produce rf pulses at 30 MHz. 
The transmitter is capable of delivering about 1.8 KW of pulsed 
power output. The receiver system consists of a pre-amplifier, 
an integrated rf amplifier and a high' gain I.E.I. amplifier. 

The overall recovery time of the receiver system is less than 
50 (J- secs. Boxcar integrator is used to sample the detected 
output from the receiver system. Electronic scanning of the 
pulses enables one, to record the recovery of magnetization 
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automatically on a strip chart recorder as a function of time. 

fhe sample holder used in ihe experiments is made of Be-Cu alloy. 

The pressure seals are satisfactory upto 3000 psi at room 

temperature. CH^, He, He, 1^2 and CO 2 gases used were of research 

grade and were obtained from Matheson Gas Co., Ar was obtained 

from Indian Oxygon Co., and was of 9S*97^ purity. 

The theory used in the analysis of the experimental 

results of this work is described in Chapter III. following 

1 Q 

the treatment given by Bloom et al expression relating T^ to 
in tram-olecular dipolar interaction and spin-rotation interaction 
have been arrived at, and considering the relative contributions 
of these two Interactions to the relaxation rate, it is shown 
that in CH^, spin- rotation interaction is the dominant mechanism 
of relaxation. 

A general procedure to evaluate the average lifetime of 

a molecule in a J state is outlined for a given anisotropic 

potential, first order perturbation theory is used to calculate 

the transition rates. The attractive part of the anisotropic 

intermolecular potential in case of the pure gases CH^, Cf^ and 

Sif^ is written in the most general form of multipole interaction 
47 

following Gray and is given by an oc to pole-octo pole interaction 
the isotropic potential being either a hard- sphere potential or 
a Bennard- Jones potential. 

In polyatomic mixtures CH^-H 2 and CH^-C02» the attractive 
part of the anisotropic potential is given by the interaction 



between the electric octopole mornent of CH^ and the quadrapole 

moment of Kg or COg. In CH^- Inert gas atom the attractive 

-7 

anisotropic potential has a term proportional to r as derived 
50 

by Buckingham for the interaction between a spherical and 
a tetrahedral molecule. The isotropic potential is considered 
to be adequately described only by a Lennard- Jones potential 
in all the mixtures. The repulsive part of the anisotropic 
potential is assumed to vary as r having the same angle 
dependence as that of the attractive part in all the systems 
studied. 

Chapter IV presents the experimental data of proton 
spin-lattice relaxation in pure CH^. The temperature dependence 
of can be fitted with a power law of the form 

where n = 1.47 + 0.03. This is in good agreement with the 
earlier work done in GH^. The data reported in and SiB^ 
by Armstrong and Tward also follow a 3/2 law and are presented. 
An analysis of the experimental data in these three systems 
yield the anisotropic potential parameters for the model inter- 
molecular potentials considered and also the scalar values 
of the octopole moments of CH^, GB^ and SiB^. The octopole 
moments of CH^ and CB^ are compared with the values obtained 
by other experiments and the results are discussed. 

Chapter V is a continuation on Experimental Results and 
Discussions. This describes proton spin-lattice relaxation in 
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GH^-1T2> CH^-COg snd CH^-Inert gas atom systems. The extrapolated 
values of ) to 100^ ilg, OOg or Inert gas atom, ^CH -X 

v/iiere X could bo either ITg or COg or inert gas atom give the 
contribution to (T^/^ ) due to CH^-X collisions only and are 


tabulaLted. The temperature dependence of (I^/p ) 


CH^-X 


can be 


fitted v'ith a power law of the form ^OF.-X*^ ^ ’ where 

n = 0*28 + 0.08 for CH.-He , 1.1 +0,1 for CH.-SIe and 1,25 + 0,09 

Ti r 0-e7i8-!i»- |tr S-£ 5 tiO-) ^ ^ 

for CH^-Arj^, An analysis of the data in CH^-lTg and OH^-COg using 
the theory described in -Chapter III, knowing the scalar values 
of the quadrupole moments of Fg and COg, gives the scalar value 
of the oc to pole moment of CH^. Ihe oc to pole moment obtained 
from CH^-Ng data is in reasonable agreement with the values 
reported from other experiments and also w'ith the value obtained 
from pure CH^ data presented in Chapter IV. Ihe octopole moment 
of CH^ obtained from CH^-COg system is, however, smaller by a 
factor of two. The anisotropic potential parameters and the 
octopole moments of CH^ obtained are tabulated and the results 
have been discussed. An analysts of the data on CH^- Inert gas 
atom yields the anisotropic potential parameters and the 
hyperpolarizability of CH^. The hyporpolarizability of CH^ 
obtained from the experimental data is comparable to the 
theoretical estimates. 

Chapter VI summarizes the significant results obtained 
in this work. 
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This thesis contains two appendices. Appendix A contains 
a collection of circuit diagrams used in the spin-echo set up. 
Appendix B contains a Fortran IV program written for the numerical 
evaluation of integrals of type 

OP 

I(p,n,n’) = j ( j dx [ g(x) x 

( dx [ g(x) ] ^ Jp^i(xy) )] 

t'j • ■ 

where g(x) = exp [ - p V^(x)] for the classical dilute gas, 
p = 1/KI and "^qCx) = Hard sphere potential or L-J 
Potential 

J i(xy) = Bessel function of order pi x and y are dimensionless 
P+2 

variables of integration. 


* * * 
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fflTROD UCTIOH 

X.l Intermoleoular Potential- s from _ W!P Data 

The study of intermolecular forces has received consi- 
derable attention as an active field of research over many 

years in chemical physics. Significant progress has been made 

1 2 

in the understanding of central forces between molecules, ’ 
such as are conventionally represented by Lennard-Jones type 
potentials. Generally much less is tenown about the anisotropic 
intermolecular potentials which depend on the molecular reorien- 
tation and shape. There are, of course, quite a few experimental 

3 

techniques such as sound absorption in gases, spectral line 

4 5-8 

broadening and nuclear magnetic resonance in gases' which 

have proved useful in providing valuable information about the 
anisotropy of intermolecular interactions. These techniques have 
the advantage that the,'.^ are sensitive to the anisotropic inter- 
molecular forces and relatively insensitive to the much investi- 
gated isotropic forces. 

Thus nuclear spin relaxation studies in gases provide 
one of the important methods to probe the details of the aniso- 
tropic intermolecular potential which causes transitions between 
the molecular rotational states. The spin-lattice relaxation 
time (T.|), which is the rate at which the nuclear spins approach 
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thertnal equilibrium with their surroundings, can be measured 
expo rimen tall;?. The nuclear spin transitions which bring the 
system to equilibrium are caused by the local fields set up at 
the nuclear sites. These fields are rendered time dependent by 
the random reorientations of the molecule caused by collisions. 
Therefore, nuclear spin-lattice relaxation times provide detailed 
information on the process of molecular re-orientation which in 
turn give information about the anisotropic interactions beteeen 
the molecules. There have been several previous experimental 

£r -1 "I 

investigations of T^ in ITg ' and in mixtures of hydrogen 

12-15 

with other gases, where an analysis of the proton T^ of 

the ortho-Hg molecules as a function of density, temperature and 

concentration of H 2 in binary gas mixtures, yielded detailed 

information on the anisotropic part of the in termolecular potential. 

A useful step for'/zard in the direction of understanding 

the in termolecular forces would be the extension of these studies 

to polyatomic molecules and their mixtures. The temperature 

dependence of the proton spin relaxation time T^ has been reported 

1 6-1Q 19-25 

in many spheric^ and symmetric top ^ molecules and also 

24 

for CH^ infinitely diluted with He. In many pure gases where 
spin rotation interaction is the dominant mechanism of relaxation 
it ?/as found that T-j/p*^ ^“5/ ^ v/here y is the density and | 

T is temperature. This result was explained by a model in which I 

the collisions are ?;eak and the isotropic part of the inter- : 

molecular potential is given by a hard sphere potential. So far | 
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no detailed analysis has been done to obtain information on 

the anisotropic intermolecular potentials in polyatomic molecules 

and their rniichares. 

This thesis covers a study of proton spin relaxation 
in pure CH^ and in mixtures of CH^-N 2 > CH^-COg and CH^- Inert 
gas atom as a function of density (14^^ <12 amagats), compo- 
sition and temperature; the atm is to extract possible information 
on the in termolecular potentials of these systems, especially 

the anisotropic potentials; the data of fluorine spin relaxation 

18 

in pure CF^ and SiF^ reported by Armstrong et.al would also 
be considered to obtain similar information on in termolecular 
potentials. 

1.2 Outline of the Frinciples of MIB; 

The principles of magnetic resonance v;ill be briefly 

discussed here and a m.ore complete treatment can be found in 
25 

Abragam. An external magnetic field along the Z axis applied 
to a nucleus with spin angular momentum I and magnetic moment 
P =: Y I produces an interaction Hamiltonian given by 

11= - jl ^ = Y 

where Y is the gyromagnetic ratio. The allowed energies are 
= - Y ® (1.2) 

?>fhere m = -I , -1+1 +I . 
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The fractional populations of these energy levels at equilibrium 
are given by 


exp (- 1 /KT) 

> _ - - - 

® L exp(-E^/ET) 

m 


(1.3) 


where K is the Boltzmann constant and T is the absolute tempera- 
ture of the gas. The equilibrium magnetization for a system of 
N weakly interacting spins is given by 



m 

i 


P 

"m 


Y -fcm 


= N I (I+1)/3KT (1.4) 

v/here the high temperature approximation y E^/KT -^4 1 is 
used. 

The approach to equilibrium of the magnetization after 
being disturbed is often describable by the phenomenological 
Bloch equations j 


m 

dt 


= y M. X H - 


M i + M 3 M - M 
X ^ , o z 


T. 


+ 


T. 


k 


(1.5) 


where Ty and T^ are longitudinal and transverse relaxation 
times. If relaxation effects are neglected, the equation of 
motion for M in a reference frame rotating v/ith angular frequency 
CflJ may be written. 
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O M 
"D t 


Y i X (H + ) 


( 1 . 6 ) 


M is stationarv in tiie rotating frane for H = if -y * 

Tho angular freauencv ^ = - Y H is called the ’Larmor 

"00 

frequency'. If H consists of a static magnetic field and 
a time varying field H^(t) where 


ILj(t) = [ i cosciJ t + 3 sin tot ] (1*7) 

represents the applied radio-freauency magnetic field, the 
equation of motion for M in the rotating frame v.fith the X-axis 
along ( t) is. 


7? 


™ M = Y M X [ ( + ^/y ) k + i] 


= YM X K 


(1*8) 


eff 


v^here 

igff= (H^ + ^/y ) k + I (1.9) 


It can he seen from equation 1*8 that M precesses about 
pith an angiJlar frequency - Y H^^^When GO = s “eff ~ ^1 ^ 
and the magnetization therefore precesses about ILj . This 
phenomenon is nuclear magnetic resonance. 

If the time varying field is applied for a time t , then 
M will pre cess at resonance about through an angle Q - y 

and the pulse width t-gy can be chosen so that Q = 180° or 90°. 
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After the application of a 90°C'Jt/2) pulse, M p3^ cesses in the 

x-y plane about F and induces an e.m.f. in the coil Y/hich is 
*' 0 

referred to as a free induction decay. If a -n; pulse is applied 
at a time t = t after the te /2 pulse, the spins tend to rephase 
in the x-y plane at a time t=2x resulting in a free induction 
signal known as a spin echo. 

\¥hen relaxation effects are considered, the solutions to 
Siq, 1 , 5 for 

H =k. E + i 2 H, cos<A>t, E^ H (1.10) 

o 1 ’ I o 

are found to he 

M^(t) = M^y(O) [cosCoot + 0)] exp (-t/T^) (1.1 1) 

M^(t) = [sin6>3 t + 0)] exp (- t/T 2 ) (1.12) 

M 2 (t) = +[1^(0) - ry exp (-t/T^) (1.13) 

where 

M 

xy ’ X y 

Hence and ^2 determined experimentally j by 

observing ■.i_,( t) and ivl (t) (orM„(t)) respectively as a function 
of time. This thesis is concerned only with measurement of T-j 
and its interpretation. 
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1.3 A7'7o.te_ on_ zation Mechanisms'. 


In most molecular fluids, may be rritten 


1 

I. 




Rg + R^ 


(1.U) 


where ^A is the relaxation rate due to in tra-molecular dipolar 
ana quadrapolar (for I > l/2 only) interactions, i.e., those ; 

interactions v/bich transform as under rotations of 

the molecule . \vhere Cl denotes the orientation of a vector fixed 
•to the molecule, is the relaxation rate due to inter-molecular 
dipolar interacti-ons and is the relaxation rate due to spin- 
rotation interactions, whose dominant term transforms as J, 
the rotational angular momentum. In molecular gases which are not 

I 

too dense, R-p is negligible. The intra molecular di'polar and 

-D i. 

i 

spin-rotation interactions are modulated in time through fluctua- 

_ 

tions in y„ (Cl (t)) and Jft) due to molecular collisions which 

^'ul ‘ [I 

cause transitions be 'tween molecular states having different 
quantum numbers, denoted by J, K and M for spherical 'top molecules. 

The mechanism by which the spins relax was first proposed 

26 ■' 27 ' 

by Schwinger and then extended by Needier and Opechowski : 

6 . . ; 

and Johnson and ^%ugh . In case of they have obtained 1 

■ ' . ■ ■ ■ ' ''' ■' i 

expressions relating spin-lattice relaxation times to the correlati 

times of in tra-molecular interactions. Bloom and Oppenheim have | 

( 

treated the dynamics of the system with Constant Acceleration i 

Approximation (CAA) and ob'tained expressions for the correlation 
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times of the intramolecular interactions in terms of the inter- 

molecular interactions. Scattering ealctilations, have heen 

28 l3 

carried out by Kinsey et.al and Poster and Hughe imerl 

1.4 _Scope of the Present Work-. 

The theory required to interpret the experimental data 

in order to obtain information about the intermolecular farces 

is presented in Chapter III, The apparatus and experimental 

techniques that were used in obtaining the data are presented 

in Chapter II. Chapter IV describes the study of proton spin 

relaxation in pure CE. ; along with pure CH. , the data on CP. 

18 4 

and SiP^ are also interpreted to obtain anisotropic inter- 
molecular potential parameters and the respective octo pole moments. 
Chapter V deals v?ith the experimental results and interpretations 
on polyatomic gas mixtures. The analysis yields potential 
parameters and molecular quantities such as octopole moment and 
hyperpolarizability of CH^. Chapter VI summarizes the various 
results and the conclusions. Appendix A contains circuit details 
of the pulsed spectrometer. Appendix B contains a Fortran IV 
program written for the numerical evaluation of integrals of 
type 

oo oo 

I(p,nsn') = ^ jdy | [ (J dx Jp_j_i(xy) x^*^” [g(x)] x 

^ g(x)]^' ] j- 

a ■ 


(1.15) 



vAere J_, i(xy) is a Bessel function of order p;g(x) = radial 

S 

distribution function given by exp (- p ¥q(x)), p = l/KI and 
¥^(x) is given by either a Lennard-Jones potential or a hard 
sphere potential. 

Gaussian quadrature technique is used for the evaluation 
of these integrals. 


* * * 



CHAPTER II 


EXESRP I EFTAL TECHfriQUES MD APPARATUS 

II. 1 Th e Method of Spin-Echoes : 

2q 

Spin-Echoes originally developed by Hahn ^ is a very 
useful technique to measure nuclear aagnetic resonance pheno- 
mena, especially relaxation times, in a simple and direct way. 

V/hen a bulk sample containing a large number of nuclear 
spin moments is placed in a static magnetic field along 

the Z axis say, a state of thermodynamic equilibrium is attained 
with a net macroscopic magnetization M along H . The 
equilibrium magnetization can be rotated from the Z axis 

through any desired angle 0 by applying a radio frequency field 
E-| perpendicular to at the Larmor frequency of the spins. 

An intense rf pulse is used to tip the magnetization tbrough 

an angle © in a time much shorter than T^ or Tg so that the 
relaxation effects_^ while H-j is on can be neglected. The angle 
© through Y>/hich the magnetization is rotated by the rf 

pulse is determined from the relation, 

e = (2.1) 

yjhere y is the gyromagnetic ratio, is the magnitude of the 

-rf field and t,. Is the duration for which the rf pulse has 
been applied. In a static field along the Z axis, the 
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magnetization precesses around the Z axis at the Larmor 

I 

freauenev characteristic of the nuclei. Consider the situation 

when a -n ; /2 pulse is applied along the positive X* axis in the 

’frame rotating' at the rf. Follovi'ing the pulse, lies entirely 

*>50 

along the Y' axis (I'ig. la)' . The magnetization preoessing 
in the XY plane after the application of a 71/2 pulse induces 
an emf in the coil which can be detected and observed and this 
is referred to as "Free Induction Decay” (DID). As transverse 
relaxation occurs, the signal decays. In a perfectly homogeneous 
field the time constant of the decay in the XY plane would be T2* 

But the applied static field will have an inhomogeneity 
The total magnetization can he thought of as the vector sum 
of the individual macroscopic magnetizations m^ arising from \ 

nuclei in different parts of the sample and hence experiencing 
slightly different values of the applied field due to field in- 

, i 

i 

homogeneity. Thus there is a range of precession frequencies j 

centered about the larmor freouenc.w which we have taken as 

■ ' ''■ ■ ■ . 

the rotation frequency of the rotating frame. Hence the macroscopic 

! 

magnetizations m^ , kno?/n as ” iso chro mats” begin to fan out, as | 
some nuclei precess faster and some slower than the frame (Fig.lh). 
Thus because of the presence of the field inhomogeneity, the mag- 
netization in the XY plane decays in a time of the order of : 

1 /YA H^). The phase differences following the % /2 pulse can be 
reversed by the application of a % pwlse at time ^ ( 4 .< or Tg) 




p 
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after the %/2 pulse (Pig. 1c). The "isochroma ts” continue 
to process in the same direction as they were before the appli- 
cation of a -jt pulse (Pig. Id). This leads to a reclustering 
of the macroscopic magnetizations (Pig* 1e)» at a time 

t = 2'c after the %/2 pulse. The continuing movement of the 
causes them again to lose phase coherence giving rise to an 
echo (Pig. If). Thus "spin-Echoes" refer to spontaneous nuclear 
induction signals which are observed to appear due to the 
constructive interference of processing macroscopic moment 
vectors after more than one rf pulse has been applied. 

II. 2 Measurement of Spin-Lattice Relaxation Time T^ s 

Measurement of spin-lattice relaxation time T^ by the 

25 2 ° 

pulse technique is well known ’ and will be described here 
briefly. 

In the measurement of spin-lattice relaxation time, the 
system was first disturbed from equilibrium and then the growth 

of the z-component of magnetization was monitored as a function 

31 . 

of time, Bloch’s equation predicts that the growth of magneti- 
zation along the z-axis would be given as 

M^( t) = + [M^( 0) - ] exp (-t/T^ ) ( 2.2) 

where is the equilibrium magnetization. Experimentally T^ 
was measured as follows? 
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lirst the application of a % pulse tipped the magneti- 
zation Mg to the negative Z direction and its growth along the 
Z axis was monitored by the application of a n /2 -x' -% sequence. 
If X vias kept constant throughout the experiment, the magnitude 
of the echo would he proportional to the value of preceeding 
the ^/2 pulse. 

The pulse sequence Tt -x - % /2 - - Tt (t was 

.■H , 

continuously varied whereas x x/as kept constant throughout the 
experiment) was repeated after every time interval t ^ 10 T^, 

so that the spins relax and attain equilibrium in betv;eon the 
sequences. A Boxcar integrator model no. GW 160, of Princeton 
Applied Research was used to sample the detected output from the 
Ii.l.L. amplifier. The use of Boxcar integrator improves the 
ratio. In order to avoid any possible integration distortion 
effects of Boxcar integrator while measuring relaxation times, 

the integration time constant RC was chosen to satisfy the 

.32 

relaxion 

RG (( a (2.3) 

s 

where T is the repetition rate, t_ is the sampling gate width, 
a is the inverse of the rate at which the pulse separation 
varied and T.^ is the spin-lattice relaxation time. The Boxcar 
was triggered by the x /2 pulse and the sampling gate of the 
Boxcar was adjusted to sample the echo height. The magnitude 
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of tlie echo and hence the recovery of ¥;as plotted on a 
Baasch. and Barab recorder Model VOM 6. The time t between the 
first % pulse and %/2 pulse was measured with a Hewlett-Packard 
524-C electronic counter and the time was printed out automatically 
by a Hewlett-Packard 526-B digital recorder-. An event marker 
on the B^schand I-amb strip chart recorder made a mark v/henever 
the time v/as printed out by the digital recorder. The at 
various times could then be identified. Typical plots of T^ 
recovery are shown in Pigs. 2a and 2b. 

A plot of log (M„-M_(t)) versus t gives (Refer e-q-. 2,2) 

o z 

at a straight line with a negative slope equal to l/T^ . 

1^*5 The Spin-Echo Apparatus ; 

i) Generals In this section the description of a complete 

pulsed nuclear resonance spectrometer is given. Standard 

2q 55 - 

techniques as described by Hahn and Clark' wore used in 
building the spectrometer. The essential fca'tures of the spectro- 
meter are a) the operating frequency is 30 MHz b) The use of 
phase coherent detection c) the construction of a sample holder 
and d) the use of Boxcar integrator to measure relaxation times 
to a better accuracy. Circuit details are presented in Appendix 
A. A block diagram of the spectrometer is shovm in Pig. 3. All 
the major items •'Hill be described briefly. 



Time interval 



Fig. 2a ~ A Typical recovery. 
( TT - T - TT / 2 - T- TT sequence) 





30 MHz 
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ii) Spin-Scr..o Spectrometers ' 

a) Timing Circuit ; 

The timing circuit consists of a collection of 
Tektronix \7ave form and pulse generator units. The Tektronix 
type 162 waTC form generator was designed, primarily as an 
interval timer and repetition rate goncra.tor. It was used in 
recurrent mode to provide the first % pulse. The saw-tooth 
from the same generator started running down at the same time 

and was fed to a modified ultra-slow Tektronix type 162 wave 

34 - ^ 

form generator.' The modified ultra-slow Tek.162 was parti- 
cularly useful in producing the required ilowly varying delay 
of the r\./2. -tc pulse sequence v/ith respect to the first % 
pulse ?jhile sampling the recovery of the magnetization. Provision 
was made so that the d.c. level from which the run down of the 
ultra-slo?/ starts could he varieds instead always from +120 volts 
to +20 volts ( Pig. A1 ) . A reset mechanism v^as incorporated in 
the modified Tektronix type 162 waveform generator to enable one 
to stop and start the run dovra at any time instead of waiting 
for the run dovi/n to be completed of its own. 

b) Pu lse S equencers 

The pulso from the modified Tek.l63 pulse-generator 
was taken to a pulse sequencer (Pig. A2) to produce tv‘?o pulses 
?;hose separation v/as controlled by the pulse width of the input 
pulse. The pulse from Tek. 162 waveform generator and the 
second pulse from the pulse sequencer were added together in a 
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mixer circuit. These two pulses were taken to the 180^ channel 
of the pulse width generator whereas the first pulse from the 
pulse sequencer went to the 90° channel of the pulse ?i;idth 
generator, 

Bjlse-k^id th Generator and P ulse- Am plifie r* 

The circuit diagram of the pulse width generator is 

given in Fig. A3. The two channels of the pulse width generator 

3R 

were identical and Y;ere designed using Phantastron circuits. 

The screen and the suppressor voltages of 5725 were obtained 
from a divider arrangement These resistors were chosen 

such that in the quiescent state the suppressor grid was suffi- 
ciently negative so that no plate current flew, all the space 
current going to the screen. Triggering was done with negative 
pulses applied to the plate and hence fed to the grid through 
the capacitor C. A negative trigger so applied, reduced the 
cathode current and consequently raised the screen voltage. Thy 
rise of screen voltage was transmitted to the suppressor to bring 
the suppressor above the point of plate current cut off. The 
trigger size was large enough to start the regenerative action 
but not so largo as to cut off the tube current. When the tube 
is driven below cut off, the grid voltage will rise initially 
with a time constant EC into the conducting region, and there 
will be a delay betv/een the application of the trigger and the 
start of the sv;eep. The pulse-width was controlled by adjusting 
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the voltage A cathode follower was interposed between the 

plate of the amplifier and the capacitor C in order to have 
shorter fall times for the pulses. The last stage was a mixer 
stage, the output of which was taken to a pulse amplifier for 
further amplification. The range of the pulse width could be 
varied from 0-35 p sec continuously. It is possible to increase 
the range of the pulse width still further by changing EC appro- 
priately. The circuit of the pulse amplifier is shown in Fig. A4. 
The amplifier stage was driven from cut off to saturation. The 
output stage ¥;as a cathode follo¥;er. A beam power tube was used so 
that the loaded pulse amplitude does not fall below 110 volts. 

d) Sadiofre qu ency Transmitter ; 

The schematic diagram of the 30 MHz transmitter is 

shown in Fig. 4. The 30MHz transmitter was designed to deliver 

high pulsed power output, and short rf rise and fall times. To 

meet the power requirements class C operation was preferred in 

certain stages follC’,vtng the rf gate. A crystal controlled 

oscillator (Fig. A5) at 10 MHz was housed in a rf leak tight 

36 

copper can. In the master oscillator design which is of 
modified Golpitts type, it is noted that the capacitance connected 
finm base to emitter is several times more than that connected 
from emitter to ground, thus effecting an impedance transformation. 

The ceramic triode 7077 was used in the grounded grid 
configuration for gating the 10 MHz rf from the source follow’er. 
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The rectangular pulses from the pulse amplifier were used to 
turn the triode 7077 on, for the duration of the pulse widths. 

The gated 10 IIHz output was taken to a 10 MHz amplifier through 
a cathode follower^ then it was tripled and amplified by a pushr- 
pull stage. The output from the push-pull stage was taken to a 
final po?/or amplifier stage at 30 MHz, The power amplifier 
produced about 1 ,8 177 of pulsed power output (Pigs. A6, A7 and 
A8). 

The output from the master oscillator went through a 
phase-shifter to a tripler (Pig. A9). A commercial delay line 
Ad— Yu Type 559® was used as a phase-shifter which prorided a 
delay of 0°-180°. The reference voltage obtained from the tripler 
output was coupled to the last stage of the high gain amplifier 
(L.E.L. Receiver) for phase sensitive detection, 

e) Sample Coil ; 

A single coil which allows optimization in the trans- 
mitting and receiving modes was used as a sample coil. It was 
made of approximately 12 to 15 turns of 20 SFG- wire and v;as about 
1^” long. The insulation on the wire was completely stripped 
off and the wire was cleaned thoroughly before use. A thin- 
v/alled glass tube was fitted tightly inside the pressure vessel 
to insulate the coil from the vessel. The sample coil v^?as wound 
to fit exactly inside the glass tube. The space outside the 
sample coil was plugged so that the entire gas sample was 
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confined to the coil. In the transmitting mode the sample coil 
?/as required to have low Q in order to avoid ringing of rf 
pulses after their application. While picking up the nuclear 
induction decay it was required to have high Q. A loose coupling 
scheme was adapted and the schematic is shown in Fig. 5. 

Fr e- Am pi if ie r ; 

The circuit diagram of the pre-amplifier is shown in 
the Fig. A10. It is similar to the design by Clark'. The 
pre-amplifier ?7as required to have high sensitivity and recover 
quickly from large overload' signals. Use of the tube 7722 at 
the input stage provided low input noise as the ^ ¥;as high. 

A tuned circuit at the input enabled maximum input to the 
pre-amplifier when the signal was picked up. The crossed diodes 
at the plate of the first stage prevented loading of the subsequent 
stage. The output stage v?as a cathode follower, with an output 
impedence of about 80 ohms. The voltage gain of the pre-amplifier 
was about 10 

g) Integrated rf Amplifier ; 

The RCA 3005 integrated rf amplifier was used following 
the pre— ampl if ier to further increase the voltage level of the 
signal before it vms fed to the high gain amplifier. The input 
and output inpedances of CA 3005 were matched for 50 ohms. The 
design of CA 3005 amplifier was taken from the RCA application 




Fig. 5 -Schematic:coupling of RF head. 
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note ICM-5022 and tiie circuit elements used are shown in 
5'ig. A1 1 . The integrated amplifier had a power gain of about 
25 db. 

h) High-Gain Recei-yer ; 

The rf receiver was a L.E.I. amplifier of the varian 
associates, model Ho, 1M2-30-1 5-50. The centre frequency was 
at 30 MHz, The 3 db band width was 10 MHz. Tr/o outputs v/ere 
provided, the one after detection has a gain of 90 db maximum 
and the other after a video stage has a gain 115 db minimum at 
0 volt bias. The quoted noise figure was 1.8 db. The input 
impedance was 50 ohms. 

A H295A diode detector was used in the h.E.L. amplifier. 
Ordinary diode detection is not linear over the entire region 
and the non-linearity becomes important ?/hile dealing with small 
signals. Phase sensitive detection v/as used ¥/hich allowed 
operation in the linear region of the diode. The diode detects 
the sum of the reference voltage and the signal. The signal 
was kept below l/lO of the reference voltage to avoid distortion. 

An experiment was performed to find the biasing voltage of the 
diode for linear detection. Tbr a given setting of the spectro- 
meter, the reference voltage, biasing the diode detector was varied 
and the output signal level was recorded each time. The EUR 
signal strength as a function of the reference voltage is 
shown in the Pig. 6, It is obvious that for linear detection. 




o cr> GO t^ CD iO ^ 

^ indino 1VN9IS 


Bias Voltage 

Fig. 6 - Diode detector response. 
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the reference voltage biasing the diode detector oust be kept 
between 1.8 and 2.8 volts negative. 

i) Sample Holder and Gas Handling System ? 

The high pressure sample holder (Pig. 7) consists 
of a thick walled Be-Cu vessel (o.d. 2.7 cm and i.d. 1 cm) and 
Be-Cu plugs A and B. The standard high pressure sealing technique 
was adapted in all pressure sealings (1AE cone fitting). The 
seats Yiexe machined to have a taper of 60® while the plugs 1 
and B have a taper of 59°. The retum lead of the sample coil 
was fastened to the plug B by means of a screw. The plug A was 
first tightened to make the seal and then the plug B was sealed 
in its position by tightening the nut E, The rf head was sealed 
against pressure at room temperature v/ith a teflon washer. The 
pressure seals were satisfactory upto about 3 >000 psi at room 
temperature. 

The pressures were measured using Bourdon pressure gauges 
which ?;e re calibrated using a dead weight gauge of Ruska Instru- 
ment Corporation, Texas. 

The gas handling system is shown in the schematic of 
Fig. 8. All high pressure connectors used in the gas handling 
system were purchased from Autoclave Engineering Inc, 

3) Heater g 

A commercial 1000 watts non— magnetic spiral, heater 




Fig. 7 - Sample holder. 
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element was wound non-inductiTely on the outside of the pressure 
vessel using a thin mica sheet for insulation, A liquid porcelain 
called '’saurelsin" v/as used to hold the windings in place, D,C, 
was used to heat the sample holder since small 50 Hz noise ?ms 
observed when A,C. was used. About 80 watts of electrical power 
had to be dissipated in order to maintain the sample holder at 
about 600°K, The pressure vessel along ?;ith the heater was kept 
inside a vacuum jacket to keep the conduction and convection 
current losses to a minimum. The leads from the heater were taken 
outside the vacuum jacket through Kovar seals. 

k) Temperat u re Measurements ; 

The temperature was measured using calibrated Chromel- 
Alumel thermo-couple embedded on the top portion of the pressure 
vessel away from the heater element. The thermo-couple leads 
were also taken outside the vacuum jacket through kovar seals. 

A Honeywell temperature controller model no. R7086 A (millivolt 
operated) was used to control the temperature. Y/henever the 
temperature v/as changed atleast 8 te 10 hours were allowed before 
taking the measurements to ensure that the sample was at equili- 
brium Y/ith the vessel. The temperature was measured with a Leeds 
and Northrup potentiometer, type 8087. The temperature measure- 
ments were accurate to at least + 2 °. 
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Ij., 4- i) Mixing of Gases? Composition Determination s 

111 the gases that were used in the experiments except 

Argon were of research grade and were obtained from Matheson Gas 

Co., Argon 'was obtained from Indian Oxygen Co,, and was of 

99*97^ purity. A mixture of the two gases 'was prepared in a 

cylinder by first introducing the low concentration component 

and then the second component to bring the total pressure to the 

required value. The total pressures ?^ere noted before and after 

the introduction of the second component. The gases were left 

together for about 46 hours after mixing to obtain a homogeneous 

mixture. The composition of the gases was estimated from the 

partial pressures and v;as checked experimentally by comparing the 

signal strength of the mixture with that of pure CH^ under 

identical conditions. The number density of CH^ was obtained from 

57 

the compressibility data by Petroleum Pesearch Institute. 

ii) Experimental Sotermination of Composition ; 

The height of the spin-echo A( 2 t ) where x is the 

58 

time befeveen 90 ° and 180° pulses is given by 

A(2. ) = A{0) exp ^ (2.4) 

in a magnetic field of constant gradient along the Z axis. In 
C q • 2 * 4 , 
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A(2 T ) = echo amplitude at time 2 


MO) 


= ecn: 


itude at time m = o 


Tg = spin-spin relaxation time 

_D = self-diffusion coefficient of the molecule 
Y = gyromagnetic ratio of the spins considered 

6 = constant gradient of the external magnetic-field 

1 /- 

At moderate densities both and - are proportional to J , 

the density of the gas in amagat units, where an amagat is 

iq ^ ^2 

2.69 X 10 ^ molecules/cc. Assuming ~= — and neglecting 

2 2^3 ^ T 

n ) — where r was about 0.5 msec, the equation 2.4 

becomes 


A{2 T ) = A(0) {-2x /T. ) 


(2.5) 


If the number of spins per molecule is denoted by , then 


A(0)^ M(0) ?;here M(0) = 


f n-r Y ^ 1(1+1 ) 


3 KT 


( 2 . 6 ) 


or A(0) = C f 
becomes , 


where C is a constant, then eq. 2.5 


A( 2 '^ ) - C exp i-2x ) 


(2.7) 


Therefore (0 y )T,/r^ ^( 0 ^ ) 


pure CHj. 


gives Itie fractional composition 


of CH^ in the mixture. Since the mixture pressures involved 
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in all the experiments were not high, the composition determined 
by the experimental method agreed within 3 to 4”^ to the value 
determined by partial pressures. In other words though the non- 
ideality was taken into account while determining the composition 
the ideal gas la?^ vms adeauate enough to determine the same. 
Typical plots of signal strength vs density of CH^ and signal 
strength vs pressure of CH^-COg mixture are shown in Pigs. 9 and 
10 . 

All densities of the mixtures were calculated using the 
ideal gas law 

^ ^ p (in atmospheres) (2.8) 

where T is the temperature in degree Kelvin and P, the pressure 
in atmospheres. 


* * * 




-Signai strength vs. density. 





CH4-C02 

C02:83°/o CH4:17°/o 



Fig. 10-Signal strengih vs, pressure 




CHAPTER III 


THEORY 

III.1 In tiaBol ec-ai-ar Inteia_ctlon_ 5 i 

The spin-lattice relaxation rate for niost molecular 
liquids and gases receives contribution from the inter mole- 
cular and intra molecular interactions. For gases which are 
not too dense, the inter molecular contribution is negligibly 
small because of the large average separation betvi'een molecules. 
The in tra molecular interactions which contribute to T'^^ in 
dilute gases are of two types, a dipole-dipole interaction 
betv/een the different spins (I = l/2 only) within a molecule 
and a spin-rotational interaction between the spins and 
rotational magnetic moment of the molecule. In case of 
the strengths of these t?;o interactions are comparable. The 
strengths of the coupling constants in H 2 for the intra molecular 
dipolar interaction and spin-rotation interaction are 34 gauss 
and 27 gauss respectively,'" But in the case of polyatomic 

iq 

molecules such as CH^, Bloom et. al ^ have established, as 
would be shown later in this chapter, that only spin-rotation 
interaction is the dominant mechanism of relaxation. 

Collisions bety/een molecules give rise to transitions 
betv/een the rotational states of the molecules. In case of 
symmetric top or spherical top molecules, these rotational 
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states are denoted by the auantum nunbers J,K and M which are 
the values of the total angular momentu^i J , the projection of 
J on the symmetry axis of the molecule and the projection of 
J on a space fixed Z axis respectively. 

As a result of the transitions betv/een different rotational 
states, the intra-mo lecular spin dependent interactions fluctuate 
with time. These fluctuations enable the nuclear spin system to 
exchange energy with the rotational and translational degrees 
of freedom of the molecules. These non-spin degrees of freedom 
constitute the lattice at a constant temperature T to the 
nuclear spins. The nuclear spin system thereby approaches equili- 
brium at temperature T at a rate which is governed by the strength 
of the intra molecular interactions and the spectrum of frequencies 
associated with the fluctuations of these interactions. 

The eigen functions | of the rotational Hamiltonian 

of a symmetric top or spherical top molecule are characterised 
by the quantum numbers J, K and M. 

Sg j J k m) = '^00^^ 1 J K m)> (3.1) 

where 

^ J(J+1) + 'I- (1 - i (3.2) 

JK 2Iq 2 

for a symmetric top molecule. and denote the t?-'o principal 


moments of enertia 
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For spherical top molecules and hence 


_0 


( 3 . 3 ) 


7 /here is independent of TC. 

i) Magnetic Bipolar Interaction; 

The Hamiltonian for the intra molecular dipolar and quadru- 

polar interactions can be written in terms of products of nuclear 

spin operators and the Y2^(il. ). Prom the general theory of 

25 

nuclear spin relaxation one obtains the relaxation rate 
due to intra molecular dipolar and quadrupolar interactions, 
in terms of Fourier transforms correlation 

functions G-2 t ) of the Y2 ) (■fc))^ 


where is the nuclear Larmor frequency and 

0 

i <^-( 1 ^))) (?- 5 ) 


where <[ ) represents an ensemble average. The quantity 

depends on the type of interaction considered. 

Bloom et.al '' have evaluated 
spherical top molecules; 


free 


for 
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For spherical top molecules, 


Sph 




free 


c:io 

1_ 'S' (2 J+1)' 

471 — 5 — 

J=0 


+2 

y . (1 + 

n=r-2 
n+J >, 0 


■JEM 


2n 

2J+1 


r rKi r 1 

Oos 1- — (2J + n+l)| (3.6) 

where n = J'— J and is the characteristic temperature 
defined by = •h^/2E . 

The Courier transform or "spectral density" J(iA3 ) of the 

free molecule correlation function is a series of delta functions 

at nK'T^ ( 2J + n+l)/h for spherical top molecules. The effect of 

collisions among molecules is to limit the lifetime of the 

molecules in their rotational states and hence to spread the 

’spectral density' associated with each 6 -function over a range 

of frequencies of the order of collision frequency. Since the 

lowest rotational frequency is often much larger than the nuclear 

Larmor frequencies at ?/hich the spectral densities are calculated, 

the contribution to J( ^ ) from 6 functions centered at 

o 

frequencies other than zero can be safely neglected. Then the 
oscillating term in eg. 3.6 is dropped and the free molecule 
correlatton functions become independent of time. The effect 
of molecular collisions is to cause G-g ^ decay to zero 
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in a time of the order of In polyatomic molecules it 

was assumed that the effect of collisions could he adequately 
described by ¥.'riting 

) 1 ( 5 . 7 ) 

*' free 

where g 2 ( 0 ) = 1 and g 2 ('i^ ) is a function ?rtiich decreases 
mono ton ically to zero as t is increased. It is seen that in 
eq. 3.6^ [ ^-^free independent of m. Hardy'^' has shown 

that under quite general conditions, w’hich should apply to 
gaseous methane 


Gi^^C-r ) = t ) exp (i m t )) (3.8) 


where j is the rotational larmor frequency of the molecule. 
Then the spectral density of ( t ) can be written as, 

jf Ui 

-toO 



~ T% ^2 ( 3 . 10 ) 

where j g 2 ( ) exp (-it^-c ) dr (3.11) 

— oO' 

The quantity f 2 must take on a value bet^/een 0 and 1. ^Tien 
the oscillating term in eq. 3.6 is dropped, using eqs. 3.6, 3.7 
and 3.8 it is shown that fg = 0.2 for spherical top molecules. 
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Using eqs. 3.9 and 3.4 one gets, 

~A^2 [^2 + ^ 32(2 c^q-2u3j)] (3.12) 


il) Spin-I?.ot 8 ,tion Interaction; 

The spin-rotation interaction arises from the fact that 
the rotating electron cloud of the molecule establishes a local 
magnetic field at the sites of the nuclei. The transitions in 
the rotational states of the molecule caused by collisions 
modulate randomly these magnetic fields, thereby causing the 
nuclear spin transitions which bring the spin system to equilibrium. 

In polyatomic molecules, the spin-rotation interaction is 
in general characterised by six independent parameters, the tiiree 
principal values of the spin-rotation tensor and the three Suler 
angles required to specify the orientation of the principal axis 
coordinate system. In spherical top molecules such as CH^ where 
the principal axis coordinate system of the spin-rotation tensor 
coincides with the molecular principal axis coordinate system, the 
number of independent parameters is reduced to only tv-’o which 
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are denoted by C| j and . 

The spin-rotation Hamiltonian of a single spin is 


given by 


19 , 41, 42 


+1 


+1 


(-1)* y - Cad- y - ^ 


k =-1 


k =-1 


(3.13) 
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where 


~ (I + i I ) 

r 2 X y ' 


~ r 2 


(j, + i y 


I 1= — (I - i I ) 
-1 f 2 ^ ^ 


3 = -i— (J - i J ) 

f 2 ^ ^ 


[ = I 

0 z 


3=3 
o z 


0a= i («ll+ 2 C,_) 


0, = 0, - C,, 


(3.14) 


The primes in eq . 3.13 indicate that the components of I and J 
are referred to a coordinate system fixed in the molecule v/ith 
the Z axis along the symmetry!' axis of the molecule. 

Transforming I* and J’ into a SfF, Eubhard"^”^ has obtained 


-h“^ K 


m I ,m 


(3.15) 


where 


(-1)° 21 ) C(112 i m ) X 


^o s-(m+ ^ ) \x 


(3.16) 


The operators D' 


m IT* 

®1-2 


oand J and hence K., ^ are time-dependent 

[X 1 J® 


because of changes of molecular state as a result of collisions 
among molecules, from the general theory of nuclear spin 


relaxation 
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£><■? 


R 


4 = 4 J ) 


— i X 

e o .dT. 


- 


(3.17) 


where G--j']( t ) is the correlation function of ), defined as. 


(3.18) 


The bracket <(_ y denotes an ensemble average of the molecules 
at the temperature T. 

The contribution to the relaxation rate T^ in spherical 
top molecules due to spin-ro tation interaction has been calculated 
by Bloom et.al using time correlation function method^" They 
have neglected the contribution to the spectral density of ^(t) 
from -A J 0 transitions which fall in the infra-red region. 

Also they have assumed that the effect of molecular collisions is 
to cause the correlation functions associated with the freely 
rotating molecule to decay to zero exponentially. In the short 
correlation time limit (W^ "^ ) rotational lines 

are broadened by an amount lar^r than the average centrifugal 
distortion splittings (— 200 MHz), the theory in which the 

centrifugal distortion effects are neglected may be used. 
Following their theory T^ may be written as, 

^ ^0 ( 5 . 19 ) 

where 0^^^ = lo^ + (4/45) ] for a spherical top moleculel- 
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Ca and are the spin- rotation interaction coupling constants 
having values 10.4 Ji-Hz and 18.5 KHz respectively^^ for CH. ; 

J 

T ^ is the correlation tinie of the spin-rotation interaction 

corresponding to the J state. 

Considering the relative contributions of the intra 

molecular dipolar interaction and the spin-rotation interaction 
1o 

Rq 5 Bloom et.al have concluded that in the case of CH^, R_^ 
contributes 8?^ of for = 3 -tg and 20^o for =t; g ^t 
1 = 100°K and 3f“ and 8^ respectively at T = bOO'^K. y ^ and 
are the time constants of the exponentially decaying reduced 
correlation functions g^Cx ) of the spin-rotation interaction 
and intra molecular dipolar interaction respectively. It may 
also be noted tha'c there are quite general arguments for 
believing that '^2 tias between the values of 1 and 3. 

This establishes that the spin-rotation interaction is the 
dominant mechanism of relaxation in CH^ at room temperature 
and above? hence the parameters used to fit the experimental 
data as a function of temperature (300°K-600°K) would be asso- 
ciated with Rq only to a good anproximation in all the cases 


studied 
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^ 'fcQPflo Secular Pot ential and Intramolecular Correlation Tiae 

From the foregoing discussion it is obvious that nuclear 

spin relaxation cannot take place in tic- absence of collisions. 

The effect of collisions is to limit the lifetime of a molecule 

in a J state. For a two level system Bloom and Oppenheim^ have 

related the in tra molecular correlation functions explicitly 

to the intermolecular anisotropic interactions. The proton 

relaxation data in pure H2 and in mixtures of with other gases 

were Interpreted to obtain information on the anisotropic inter- 

5 8 

molecular potentials.’ In polyatomic molecules such as CH^, 
because of its large moment of enertia, appreciable number of 
rotational levels are populated even at room temperature. Hence 
the results of 112 are not applicable to polyatomic molecules. 

In fact it does not seem that the two level results of H2 and its 
isotopic modifications are applicable to any other gaseous system 
as the number of rotational levels involved is quite large. 

In the absence of a detailed molecular theory for polyatomic 
molecules, the theory of Bloom and Oppenheim has been used in a 
modified w'ay to obtain information about the intermolecular 
potential. Since there is no simple theory to relate the 
correlation time of the intra-molecular interactions to the 
in te r molecular potential, it is assumed, a priori, that t -j , 
the correlation time of the spin— rotation interaction is equal 
to the average lifetime of the molecule in that J state. The 
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average lifetime of the molecule in a J state may then be 
evaluated for a given anisotropic potential using first order 
perturbation theory and the results of Bloom and Oppenheim for 
in termolo cular correlation functions. She liustif ication for 
this assumption must then come, a posteriori and will be discussed 
later. Under this approximation, it is then possible to obtain 
the anisotropic potential parameters using eq. 3. 19* 

III. 3 G- eneral Procedure for the Evaluation of the Average Lifetime 
of a Molecule in a J State s 

Defining A(J*M'; JM) as the transition rate for the transitions 

J 11 

Jli — > J*M’, can be expressed as 

+J 

1 _ 1 
7 J “ 2J+1 

1 M=-J 
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It is assumed that the collisions are weak in the sense that 
first order perturbation theory may be used to evaluate the 
transition rates JM) . Then the transition rate A(J’M*5JM) 

is given by- 

AC J»M' I JM) = ^ y t'pjj z: A( J’K’M* I JM) (3.21) 

and 


I 


+J' 

■ ^ A(J’M’;JM) 

M'=-J* 


( 3 . 20 ) 


M’4i for J=J’ 
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-f 

A(J’K'H'; JEM) = j e"^^JJ»^ x 

J 

<^imfil 7^(0) f J'K'M'} (j‘K‘M»| Vj.(t) | Jiai)^ dt 

( 3 . 22 ) 

v^here | JKM^ are the eigen states of the unperturbed rotational 
Hamiltonian (eq, 3.1) and 

^^JJ» = % - ®j« (3.23) 

is the energy difference betv/een the rotational states J and J’ 
<(^(Jiai I '^j^(O) i J'H'M*) (J’K*M' j ’^j^(t) I JEM)^ is the correlation 

function of the matrix element of 7j^ between states J'K'M’ and 
JEM where the average is performed with the density matrix 

e“ ^ ^ /Tr [ e" P ^ ] and 

V^(t) = 7 ^( 0 ) (3.24) 

where 

is that part of the lattice Hamiltonian which contains the 
kinetic energy of the centres of mass of the molecules and 
isotropic interm olecular potential. 

i) Inteimolecular Potentials: 

The into rmolecular potential V beteeen tv/o molecules 
separated by distance r is given by 
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V = v^(r) + V^(r, ©') 


( 3 . 26 ) 


where the isotropic part V^(r) depends only on r and the 
an iso tropic part ¥^(r,Q’) depends both on r and ©’ where ©* 
denotes the orientation of the oolecule relative to a space 
fixed frame. 

Tv/o models of isotropic interaolecular potentials will be 
considered^ 

Model 1s Hard sphere potential 


Y^{x) - oc X < 1 


7^(x) = 


(3.27) 


0 X 1 

where x = r/a is a dimensionless quantity and a is a potential 
parameter. 


Model 2 % Leonard- Jones potential: 

= 4 s (x“^^- x“^) (3.28) 

^ , and a = r/x are appropriate energy and distance parameters; 


ii) Anisotropic Intermolecular potentials; 

The anisotropic intermolecular potential is in general 
written by considering the electrostatic interaction appropriate 
to the colliding pair. Specific forms of the anisotropic 
po ten tials appropriate, for the different systems will be 
discussed in the following; 
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The most general form of tlie attractire anisotropic 
interaction is given 


A% Cj 
h pl+1 


7. c( 


i.V; Y*jn) 


m:} ,m2 


(3.29) 


The orientations of O. ^ , (^2 ^ ©it relative to a 

SFP (Fig. 1l). Ihe vector connecting 0^ to Og is denoted by 
R = (R, n.) representing the intermolecular axis where d =('Q,<^ ) 
give the orientation of R, The molecules have been considered 
as having rigid non-overlapping charge distributions centred 
around origins 0^ and O2. 

A = 1^12 ’ ^ = 1^+12, 1^ and Ig indicate the order of the 

electric moment; for dipole 1-|=1 , quadrupole l-j=2 etc., 

Ai^^2 4 IX (21 + 1)1 

C = C, . = irJ_ [ ] (3.30) 

i^ ,J.2 (21+1) (21^+1) I (212+2)1 


The multipole moments Q, _ are defined by 

i- f 


Ter Y^Jr) 


3.31) 


the sum extending over all the charges in the distribution with 
r = (r,r), the position of the charge e relative to the origin 
of the distribution. For a continuous distribution of charge. 
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tlie sum in definition 3.31 is replaced by an integral over the 

charge density ^ (r). P/hen the intermolecular potential is 

considered for definite electric moments, 0, , for a fixed 1, 

'l,m ’ 

would then form the components of an irreducible tensor of 

47 

rank 1, It is shown ' for systems having axial symmetric charge 
distribution such as ^102” 

= (5-52) 

where is a constant; Cl denotes the orientation of the 

symmetry axis of the charge distribution. In case of CH^ where 

there is a symmetry the oc to pole moment is given^"^ by 


Q 


3 ,m 


= i ( 


7 


4TT 




D\(n 

m.y2 


) - 


m ,-2 


(Clj| 


(5.53) 


where 

= {G, <p .%) denotes the orientation of body fixed 
frame relative to a SIF. D’s are the corresponding rotation 
matrices^® fl' is the scalar value of the octopole moment. 

a) Spherical Top Molecules (CH^, CF^ and SiP^) s 

The attractive part of the anisotropic intermolecular 
potential in spherical top molecules such as CH^(also CF^ and 
SiF.) is given by the electrostatic interaction betv^een the 

peroiaoBnij . .octo poXqs rssidGU^’b . oo, , ■‘feh© iboIgcuIgs wlisrsas XIig 
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repulsive pci.rt is assumed to have an r dependence and the 
same angle dependence as that of the attractive part. Then 

using eqs. 3.29 5 and 3.30, the anisotropic potential is given 
hv 


¥^(r,0‘) 




hC ax) 


0(336; x 





a-j is the measure of the strength at r=a of the repulsive part 
in the anisotropic potential. 


b) and CH^-COg Mixtures; 

The attractive part of the anisot3?opic potential in 
these systems is considered to be given by the interaction 
betv;een the octopole moment of CH^ and quadrupole moment of 
^2 03? COg. 

' ' ' ' *1 ' 2 ' 

Using eqs. 3.29, 3.32 and 3.33 and assuming an r“ 
dependence for the repulsive part of the anisotropic potential 
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having the same angle dependence as that of the attractive 


par 5 


V^(r,«<) = 4it 1(S||)5 Max)T 0(325; n,B,m) v 


(3.38) 

where h(ax) = (a^/x"'^ - a^/x^) (3.39) 






= (^# /a^) 


( 3 . 39 ) 

(3.40) 


where ^ is the scalar value of the guadrupele moment of ^2 or 


CO 2 . 


iii) CH^- Inert Gas Atom;^^ 

50 

Buckingham has derived an expression for the long 
range attractive part of the tnteraolecular potential between 
a tetrahedral molecule and a spherically symmetric atom by 
expanding the Coulomb interaction in terms of multipole moments, 
the leading terms of which are. 


3U-1IJ. 

= - nvjiVj, 


r a ^ 8 a A 

I ■ ... -g-v . + Cos Cos Q Cos 0^ 

~ r r y 


( 3 . 41 ) 


oclj and 0 ^ are the polarizabilities of the inert gas atom 
and methane molecule respectively. The ane'les 0^, 0^^ and 
are as shown in Pig. 12, A is the hyperpolarizability of the 




Spherically symmctrsc 
molecule. 


I\ synr.iT.elr 


12 -Orientation of CH 4 
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me ciiaiie molecule describing the distortion of the molecule due 
to an external electric field and field gradient. and 
have the units of energy and it was suggested by Buckingham 
that they might be taken as the first ionisation potentials of the 
two molecules. 

The first term in eq,5.41 is the isotropic part of the 
long-range interaction and hence must be identical with the 
attractive part in eq.5.28. Therefore (XJ^Ug/TJ^+TJg) has been 
chosen such that 



(3.42) 


where e and a are taken from viscosity data. It is worth 
pointing out hero that this choice gives for (U^U 2 /U^+U 2 )^ a 
value approximately twice as that would be obtained following 
Buckingham’s suggestion. This difference is due to the fact 
that r"^ term is only the leading part of the potential derived 
by Buckingham while in the (12-6) potential it is the only term. 

It is also assumed that an r dependence on r for the 
repulsive part of the anisotropic potential and the same angle 
dependence as that of the attractive part. The intermolecular 
potential may then be written as, 

V = 4. e (x*"”'^— 3 ;""^) + b(ax) Cos Cos Cos (3.43) 
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where 

b(az) = x = r/a (3.44) 

and a 2 = 52 Az / a (3.45) 

In order to calculate the transition rates A(J'M‘, Jll) , 
the anisotropic potential Vj^(r,©) is rewritten in the molecule 
fixed frame (MFI’) in terms of the spherical harmonics as, 

V^(r,«) = (^)* (^3,2- (5.46) 


Using Eqs. 4,22 and 4.28 a of Eose^® the above equation can be 
transformed to a space fixed frame (SEE) as 


V_^(r,e') 



M 

(3.47) 


where ( a^y ) s.tq the Euler angles of rotation, D's are the 
corresponding rotation matrices and 0* is the orientation of r 
in the SEE. 

The details of the calculations for the different 
systems studied in order to obtain the appropriate anisotropic 
potential parameters and the related molecular quantities such 





5S 


as the scalar ^ralue of the octopole Tioaent of CH. and the 

4 

h.yperpolarizaoil Ity of frofn the experimental data are 
presented in Chapters IV and V. 


* * * 
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RESULTS MD MUJjYSIS (pure &ASES) 

Results ;(CH_^. CR^ and SiR^) 

Proton spin-lattice relaxation tiues,!^ I'^ere measured 
in pure CH^ as a function of density ( 1 f ^ 12 amag^ats) in the 
temperature region 300— 600°K, At each temperature 1^/^ was 
plotted as a function of ^ and the extrapolated value to j’ = 0 
was taken as true value. In general it was observed that 
is proportional to density in the region 1 < |X 10 amagats, but 
slight deviations occur for densities greater than 10 amagats, 
approximatoly . Since at the moment it is believed that these 
deviations are due to surface effects, all the measurements were 
restricted to about 10-12 amagats only. The value of 
obtained at 300°K (21.8 msec/amagat) is in good agreement with 

51 

the value obtained from the lo?/-density data of Beckmann et al. 
They have reported that it is likely thot the earlier results of 
Lalita and Bloom^^ in pure CH^ were influenced by a small contri- 
bution (— 10?fO to the relaxation rate due to surface effects. 
Gerritsma et. al^^ have reported T^ measurements in the region 
108-549 amagats. They found that T^ is proportional to density 
upto about 400 amagats and obtained a value of T^/y? 22 msec/ 

amagat at 298°K. One important difference betvfcen the pressure 
vessel used by Gerritsma et.al and the vessels used by others 
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is that G-criitsrna et.al encased tbeir rf coil in glass whereas 
others had their rf coils Mmersed in the gas sa-nple. Also 
Gerritsma et.al wont through an elaborate cleaning procedure 
for tiieir sample containers. If the surface effects are coming 
from the collision of the gas molecules v^ith the surface of the 
vessel, made of glass in all the cases, Gerritsma’s results 
should have been affected to the same extent as those of others. 
Since they are not, it is likely that the surface effects might 
be coming from the collision of the gas molecules ?/ith the coil. 
Typical plots of T-| versus ^ and corresponding ^-j/^ vs f axe 
shown in Pig, 13. The temperature dependence of ^rure CH 

is shown in Pig. 14-. The data can he fitted ¥^ith a power law 
of the form T^/^:>'^ T~^ ?;here the T is the temperature. The value 

of the Index n is found to he 1.47 + 0,03. 

18 IQ 

Armstrong et.al have studied -^P relaxation in CP^ 
and SiP^ in a limited temperature region 230-370°K. The tempera- 
ture dependences of T^/y in these systems are found to be the 
same as in CH^ with the index n taking the value 1.5. The present 
values of T^/y at different temperatures in CH^ along with the 
values of (Ty/O ) in CP^ and SiP^ reported by Armstrong et.al are 
tabulated in Table 1. The effective cross section ^ for 

spin-rotation intoraction along with the corresponding kinetic 
cross sections ratios ^ ‘^eff ^^4* ^^4 

and SiP^ are tabulated in Table; 2. 



Fig. 13 - Typica 
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TABLE 1 

"VALUES OF CH^jCE^ AFI) SiE^ 


Sys tern 




SiE^ 


Temperature °K (T^/^ ) msec/aoagat 


500.0 

21.80 

400.0 

13.80 

500.0 

10.20 

600.0 

7.80 

244.5 

2.65 

300.0 

2.03 

534.0 

1.59 

571.5 

1.36 

251.5 

26.00 

295.0 

18.60 

520.5 

15.50 

537.5 

13.60 
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TABLE 2 


COM PARISOI-T OE KBTETIC CROSS SECTIONS A>T D COLLISION CROSS 

SECTIONS 


System 


Kinetic 
section 
a _ 
kin 


cross 




Temperature, 

°K 


Collision cross 
sections 
cr (\ 2 ) 
eff ■ ^ 


Mean Number 
of collisions 


( 


a / A 

kin'^ effy 


45.27 

CH^ 


CP^ 69.5 


SiP^ 88.0 


300.0 

19.46 

2.32 

400.0 

14.27 

3.17 

500.0 

11.80 

3.84 

600.0 

9.88 

4.58 

244.5 

40.5 

1 .72 

300.0 

34.3 

2.03 

334.0 

28.4 

2.45 

371.5 

25.6 

2.71 

231.5 

86.5 

1 .02 

295.0 

73.1 

1.20 

320.5 

61.2 

1 .44 

337.5 

54.5 

1.61 
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IV. 2 Interpretation ; 

As it was pointed out in the introduction, the temperature 

dependence of nuclear spin-lattice relaxation times, , available 

until now in pure gases for which the snin-rotation interaction 

is the dominant interaction, can be fitted by a power law 
-3/2 . ’ “ 1 J 

T wrhere I is the temperature and is the density. These 

results for such spherical and symmetric "top molecules are in 

accord w?ith the theory based on hard— sphere intermolecular 

potential sind wreak collision approximation; In fact such a 

temperature dependence results irrespective of the form of the 

anisotropic potential assumed. Hence such data were used so 

far only to establish the dominance of spin-rotation interaction 

and the adequacy of the hard sphere potential to describe the 
16-23 

isotropic part and no attempt has been made to obtein the 

parameters of the anisotropic part of the intermolecular potential. 

In the following analysis the temperature dependences 

1 R 

of proton in CH^ and fluorine in and SiP^ have been 
interpreted to obtain anisotropic potential parameters and 
corresponding scalar value of the octopole moment. The isotropic 
part of the intermolecular potential is considered to be described 
by a hard-sphere potential as given by ea.3.27. Since in this 
case a T“^^^ dependence results independent of the form of the 
anisotropic potential, it is not possible, knowing only the 
temperature dependence of --j/p > determine uhe strengths of 
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bo til tbo* repulsive and attractive parts of the anisotropic 
potenbial. Hence only an attractive part of the anisotropic 
potential given in eq. 3.34 will be assumed. 

In viev^ of the fact that 12-6 Leonard- Jones potential 
is more realistic than a hard sphere potential, the data on pure 
gases GH^, and SiH^ vrould also be analysed using a Leonard- 

Jones isotropic potential for the isotropic part. The form of 
Lenuard-Jones potential is given in ea. 3.28. In this model, 
in order to explain the observed temperature dependence of 
in pure gases, it is necessary to assume as it is explained in 
the following discussion that the anisotropic potential contains 


both the repulsive and attractive parts. 

The effect of the anisotropic potential (eq. 3.34) is 


to produce transitions of molecule 1 (eg. CH^) from the state J* , 
K’, M’ to the state JUl while the molecule 2 (eg. another CH^) 


s irau 1 tan e 0 u si y u n ci e r go e j 


transition from 


to J' 


K" , M". Denoting the transition rate for this process by. 


1((JM, J’K’M’) (J''K"5'["«, J"* K'", 


M'" ) ) one could write 



(4.1) 
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Tile tra.nsition rates AiJW* -k^ i 

AvojMd, d m j to be evaluated are those 

associated vvnth molecule 1. 

Using the anisotropic potential as given in eq.3,34 


and using eqs. 4.1 and 3.22 we write. 


A(JKBI; J’K'M') = Z > T < T ( jTt X 33264 x .r x 


M" M'" M' 
K” K'" 


JJ’ j J”J'” 


[c{336;II^M2T«)] 


■J" 


» > 




2! 






)) j j J’E'H')) 


"4- 


((jEi 1 d|* „( n 1 ) 


where 

j(U) 


"1 
4 


J’K'M'/ 


1 


(4.2) 


1 


JJ';J"J"'^ -^2 j 


X J 

e JJ';J''J’'' 




*c*0 


<^b(s(0) )Y6^_^( A (0) )bf 1( t)Yg^(A(t) I 
(4.3) 

The bar denotes an average over the equilibrium ensemble and 


^ ^JJ’lJ"J"’ " (Ej-Ejr ) + 


(4.4) 


and 


h'" " xrT2FJire-^"'‘^‘''+’'®o''’^ 

J'" 


(4.5) 
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where is the characteristic temperature of the molecule in 
degree Kelvin. In Eq. 4.2 C is a clebsch-ffordon coefficient. 
We use normalized symmetric top wave functions 



(2J+1)2 


'{ 8 % 



(4.6) 


to evaluate the matrix elements in eq. 4.2. After carrying out 
the summations over ^ 2 , 1 " and M'" using standard methods of 
Rose and using eq. 3.21 one obtains 


A(MiJ'M') 


TCX3 3264 X 52 n p . 2J"+1 v , 2J+1 s 3 . 
Z- ^325x49x49 ^ ^ 2 J ’+1 

C(J3J' ^ 3 ( 10 ) (4.7) 


The spacing between the rotational levels of CH^ (also CP^ and 
SiE^) is such that at temperatures of our experiment we can 
assume 3 (W ) 3(0)^ Row using the expression for 1/ “Cji 

ea . 3.20 find summing over J" with J" = J"’ + 5 "to J‘”'-3 we 
arrive at 

1/ -Tj' = 1337 3(0) (^*8) 

Ihe function 3(0) has been evaluated by Bloom et.al using 
constant Acceleration Approximation (CAA)? It is given by 


_fa^(2Trp|i )^^^ I(p) 



3(0) = 


(4.9) 
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where 

I(P) = -i-(p,n') + I(p,n) - 2a^a2l(p,n,n’ ) 

(4.10) 

and 

f ■' r ^ 

I(p.n) = jay!] ax [g(x)]^ 

0 “ (4.11) 

oO oo 

I(p,n,n') = j ‘3y| jf j dx[g(x) ] Jp^i(xy)x^'^^”” )x 

( j dx [(g(x) Jp^i(xy) )j 

^ ( 4 . 12 ) 

In above expressions p and a are the reduced mass and distance 
parameters of the system under consideration resrectively . ^ is 
the number density of the mixture, p = l/KT where K is the 
Boltzmann's constant, T, the absolute temperature^ x=r/a and y 
are the dimensionless variables of integration. Also 

g(x) = exp |^-■VQ(x)/KT J (4.15) 

is the radial distribution function and J._,i,(xy) is a Bessel 

P+ P 

function of order p. Bor octopole-octopole interaction p 
takes the value of 6; the radial dependence of the attractive 
part of the anisotropic potential is given by the index n = 

The repulsive radial dependence n' is assumed to be 12. 

Using eqs. 3.19, 4.8, and 4.9, can be written 

asj ■ ■ 




70 


T. 


4.11:' 


'eff 


1337.0 a^f 2 TEpU 


1(6) 

(4.14) 


The values of 9^^, 0^^^, a and \i are tabulated in Table 3 for 
the three rnolecules CH^, and SiP^i f^=2.69 ^ lo'^^ mols/cc 
and 1(6) is given by eq . 4 . 10 with p=6, n=7 and n’=12. 


Since in case of the pure gases CH^, CP^ and SiP^ the 
experimental data fits a 3/2 law it can be seen from oq. 4.14 tha.t 
^ ^(^^^ ezpt terfiperature independent. The temperature dependences 
of l(p,n) and I(p,n,n’) in eqs. 4.11 and 4.12 come only through 
the radial distribution as given in eg. 4.13^ If ■7 q(x) is given 
by a hard sphere potential I(p,n) and I(p,nin*) would be indepen- 
dent of temperature. Therefore it is not possible in this case 
to deterDiine both the parameters and 02 in eq. 4.10 from only 

the temperature dependence of T^/^2 as can be seen from eq. 4.14. 
Hence it was assumed that the anisotropic potential in this model 
is purely attractive. The anisotropic potential parameter a 2 
was obtained by fitting the expression given by eq. 4.14 to the 
experimental values of ^-j/^ each system over the entire tempera- 
ture region using least squares criterion. Hor the hard-sphere 
isotropic potential the value of 1(6,7) is tabulated by Bloom 
et.al and is 0.0104. The anisotropic potential parameter a.^ 
for the three systems along ?i’ith the scalar values of the 
octopole moments obtained using the equation ^ 
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TABLE ^ 

VALFES OE a MB 


System 

0°(K) 

^eff 

(MHz) 

a 

( 10~®cm) 

(gms) 

OH^ 

7.680 

158.60 

3.80 

13.40 

GE^ 

0,270 

39.80 

4.70 

73.60 

SiE^ 

0.1970 

5.94 

5.59 

87.00 
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-^1 - (a^ (4.15) 

for the hard sphere isotropic potential are tabulated in 
Table 5. 

If the iso cropic part is given by a Leonard- Jones 
potential I(p,n )5 and I(p,n,n’) would be temperature dependent. 
Though I(p,n) and I(p,n,n') are temperature dependent it is 
possible to make 1(6) temperature independent by choosing the 
strengths of the repulsive and attractive parts of the aniso- 
tropic potential appropriately and explain the observed temperature 
dependence in pure gases. Por the isotropic Lennard-Jones 
potential, the values of 1(6,12), 1(6,7) and 1(6,7,12) were 
evaluated as a function of temperature given by pe and are 
tabulated in Table 4. The anisotropic potential parameters a^ 
and a .2 for this realistic model obtained from eqs. 4.14, 4.10 and 
Table 4 ^ using a least squares fit of the experimental data over 
the entire temperature region along with the respective octopole 
moment obtained using eg,. 4.15 are tabulated in Table 5. Ibr 
purposes of comparison Table 5 also contains values of octopole 
moments obtained from other techniques. 

IV.5 Remarks on the Humerical Integrations ; 

The integrals given by eqs. 4.8 and 4.9 fo^ s given 
p,n and n* were evaluated numerically by Gaussian quadrature 
techniques on an IM 7044 computer. The subroutines were 
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TAEI.E 4- 


, Kp ^n) mi) K^.n^n*) EOR THE CLASSICil, 
DILUTE GAS FOR LElillLARI)- JQffES ISOTBOPIC HDTEImTIAI 


|3e 

1(6,12) 

1(6,7) 

1(6,7,12) 

0.621 

0.01677 

0.02571 

0.01958 

0.505 

0.01861 

0.02561 

0.02051 

0.480 

0.01917 

0.02566 

0.02084 

0.360 

0.02346 

0.02669 

0.02340 

0.288 

0.02827 

0.02818 

0.02630 

0.240 

0.03350 

0.02986 

0.02950 
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developed for the Gaussian quadrature and for the generation 
of the half-order Bessel functions. Convergence and accuracy 
of the integrations were checked by extending the upper limits 
of the variables of the integration and increasing the i»int 
densities until results of successive integrations agreed to 
the desired accuracy. Appendix B contains the details. The 
accuracy of the generated half-order Bessel functions was tested 
by a comparison with standard tables. Also the results are in 
agreement for the selected integrals I(p,n) and I(p,n,n') 
tabulated as a function of temperature by Bloom et al®. It 
is to be noted that Bloom et.al have tabulated the numerical 
values of these integrals for the case I(p,n) for lennard-Jones 
isotropic potential where p = 0,2 and 4 and for I(p,n,n') where 
p=2 only. 

IV, 4 Discussion s 

Ti/vo model intermolecular isotropic potentials, a hard- 
sphere potential and (12-6)1-1 potential, with appropriate 
anisotropic potentials have been considered for the analysis 
of the experimental data on pure gases CH^,CP^ and SiP^ and the 
results are tabulated in Table 5, Both these models successfully 
explain the temperature dependence of the spin relaxation data 
in the region of interest. It is seen from the Table 5 that 
the ootopole moments of CE^ and obtained from the hard sphere 
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.odel are so.ewhat .lese' (2^30^) ^ ^ 

.odel ecnaistlng of a lennard-Jones potential and an anlsotrorlo 
potential haying both the attractive and repulsive parts. The 
scalar values of the octopole mooients of OH^ and CP^ obtained 
froa the latter »del are in better agreement with the reported 
values from other techniques. However in case of SIP^, a hard 
sphere model predicts a value of 5.3 x 10 " esu ctn^ for the 
scalar value of the octopole moment, whereas the hennard-Jones 
potential rno del gives a value of 4.14 x There is little 

to choose between these two models as no experimental value of 
octopole moment is known from other techniaues. 



* * * 
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V.2, CE^-Hg and CH^-C 02 Mixtures: Results; 

Plo ts of P-]/^ versus are shown in Pigs. 15 and 16 
for CH^-Eg and CH^-GOg systems respectively. Figures 17 and 18 
represent T^/f as a function of Ng 

different temperatures. Solid lines in all the plots indicate 

least square fit of the data. The error bars in Figs. 17 and 18 

are the standard deviations in ). The extrapolated values 

of (1-^/f ^CH -N ) 0 jj _QQ 3.re tabulated in Table 6 . 

The temperature dependence of the extrapolated values of (T^/^ ) 

for the two systems are shown in Figs. 19 and 20. The data 

can be fitted with a power law of the type (l-j/? )qjj -TT” 

4 " 

where n=0.87 +0.14 for CH^-Ng and n=0.91 +0.1 for CH^-COg. 

The effective cross sections spin- rotation interactions 

along with the kinetic cross sections are presented in Table 7. 

The ratio ( ^ *^eff^ give the mean numler of collisions 

required for randomizing rotational angular momentum. 


Y. 3 Interpretation ; 

The anisotropic intermolecular potential as given by 
eq. 3.38 is used to calculate the average lifetime of a molecule 
in a J state. Using eqs. 3.22, and 4.1 the transition rate 
could be written as, 




CH 4 56% CO 2 44% 



CH 4 17% CO 2 



SL. 

O 


398-5 K 




-9 


630°K 


— ^ -8 12 

Density in amagats 

Fig. 16 - Plo ts of Til? versus f . 



Fig. 17 - Dependence of T]/f on N 2 concen tro I ion a t different 

tempera tu res. 




Percentage of CO 2 

F'ig.1 8 - Dependence of 1)1^ on CO 2 concentration at different 

temperatures. 



Fig. 19 -Dependence of (Ti/j temperature. 
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TABLE 6 


EXTRAPOLATED 

YilLlES OR (T^/^ ) TO 100^ N 

2 illD CO 2 at various 

TSiEERATDRES 

Temperature 

Standard Deviation 

^CH^-COg - 
Standard Deviation 

300 

13.94 + 1.21 

21.97 + 1.35 

400 

11.15 + 0.84 

17.14 + 0.96 

500 

9.00 + 0.67 

13.86 + 0.74 

600 

7.63 +0.58 

11.70 + 0.61 



o 
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A((JSiI| J'K’M* )C J”'M” • ) ) = 21 2, g52xl6 v . 

^ 55 J' J”J ” » ^ ^ 


^1®2 


[C(325; x 

J'K'M')>|2 +j((jKM| J'K'JJ ')>j ^ 


(5.1) 


where 


D# 

1 ^ 


J''») = “2 i < 3 Te 

^ J 




<’ b(r(0)) y| (il(o)) b(r(t)) Y. ^(n(t))\ 

^ ^ D fin / 


Using Rose^® 


(5.2) 


[C(J‘’»2J";00) ] -^ 

(M'*’ + = M") (5,3) 

f<^(jKI«l ^2(^-1) I J'K’M’))P = (|fT^^)[c (J 3 J’?M,-M^,M») x 

[C(J3J‘| K,-2, KO ] ^ 

(5.4) 

It is to be noted that the sum 

[c( J" ’ 21 "; 00 ) f = 

it 

u 


1 


(5.5) 
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Substituting eqs. 5.3, 5.4 in eq. 5.1 and using eqs- 5 . 2 , 5.5 
and carrying out the sumaations using standard, methods, it i 
written 

A(J1,I;J'M') = (im-JL) 3(o). JJ-H ^ [0(J3J',M,-o.^)] ^ 


(2J’+1)' 


(5.6) 

with the approximation 3 (t^ ) crT 3 ( 0 ). How using eq. 3.20 
the average lifetime of the molecule in a given J’ state can 
be written as 


245 

864 X 16ir 3(0) 


(5.7) 


where 3 ( 0 ) is given by 4.9 using eqs. 3 . 19 , 5 . 7 , (T^^ )qjj 

4 

at a temperature T is given by 

216 X I6ccf ^ a^(2 )^' 






I(P) 


(5.8) 

where. X could be I., or COg and a = 7.68°K/[D°K; Ifp) is given 
by equation 4 . 10 with p=5, n =6 and n’=12. The values of the 
integrals 1 ( 5 , 6 ), 1 ( 5 , 12 ) and 1 ( 5 , 6 , 12 ) at different temperatures 
are tabulated in Table 8 . 

Using eq. 4.10 and the numerical values of Table 8 the 
parameters a^ and a 2 were obtained by fitting the expression 
given by eq. 5.8 to the experimental data of (7^/(7 ) over the 
entire temperature region using least squares criterion. The 



GAS MIXTURE FOR LBTITARD-JQHES ISOTBOPIC POTSTiijas 


pe 

1 ( 5 , 12 ) 

1 ( 5 , 6 ) 

1 ( 5 , 6 , 12 ) 

0.19850 

0.04536 

0.03926 

0.03760 

0 . 2 U 20 

0.04235 

0.03844 

0.03580 

0.26775 

0.03370 

0.03646 

0.03142 

0.29200 

0.03119 

0.03585 

0.02999 

0,35040 

0.02676 

0.03491 

0.02749 

0.35700 

0.02637 

0.03484 

0.02727 

0.43800 

0.02248 

0.03436 

0.02520 

0.58600 

0.01919 

0.03413 

0.02350 


0,02350 
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octopolo niomant is tben given by 




i L 



^GH4-X 


(5.9) 


v/liere ^ ^ is the scalar value of the guadrupole uonient of 
IJg ox COg. 

The anisotropic potential parameters a^ and ag, for the 
systems and CH^-COg are tabulated in Table 9 along with 

the scalar value of the octopole moment obtained. 


"V . 4- Discussion of the Results s 

It can be seen from the Table 9 that -O. oblained. from 

the analysis of CH^-Hg system is in reasonable agreement with 

the reported values. The value obtained from CH^-COg systetr 

is, however, small by a factor of 2 approximately, Btllingslay 
54 - 

and Krauss have pointed out that the hyperpolarizability of 

the molecule is not taken into account while obtaining the 

quadrupole moments of the molecule and hence the recotamended 

values of the quadrupole moments are over estimated by 20^ 

approximately. If the quadrupole moments of and CO„ obtained 

55,56 

from birefriengence method are used after correcting for the 
hyperpolarizability by decreasing the quoted values by 20 ^ 0 , 

_X/1 

we obtain = 2.64 x 10 esu cm'^ from CH^-Hg and 
XI =1.33 X lO"^'^ esu.cm'^ from GH^-COg. This value of 
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2.64 X 10 is in good agreement with the -value obtained from 
the analysis if 2.^ data in pure CH^ presented in Chapter IV, 

Analysis of the data for radial dependences such as 
-15 -8 

r and r for the repulsive part in anistropic potentials 
yield values of ochopole moments of CH^ 30-40^’ smaller than 
the values obtained from an r~^^ dependences for toth CE^-lg 
and CH^-C02 systems. This suggests that collision cross-sections 
are sensitive -to form of the repulsive potential. In the 
analysis it is assumed that the angle dependence of the repulsive 
part is the same as that of the attractive part without any 
justification. It might be that they are sensitive to the 
angle dependence also. COg being a triatomic linear molecule 
longer in size than it might be more sensitive to the angle 
dependence than II 2 molecule. 


V.5 GH^- Inert Gas Atom; Results; 

Typical plots of T^ versus ^ are shown in fig. 21 
for 77.3/^ CH^ and 22.7^ Ar mixture at three different tempera- 
tures. figs. 22 , 23 and 24 represent as a function of 

inert gas atom concentration for the three systems (He, He and 
Ar) at different tempera-fcures. The extrapolated values of 
(T^/^ ) to 100^ He, He and Ar are tabulated in Table 10. The 

^ where X could be either 
4 "^ 

He, He or Ar is shown in fig. 25 . In all the cases the data 


temperature dependence of 





00 o (X> ' o 

OJ 04 CN — 
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Fig. 24 ~ Depend^n 
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TABLE 10 

V/iLUES (0? ^CH -X i STMDAHD DEVIATION AT 

~"’4 

DIPEEBKTT TEMEERATOHES 



298 7.7 + 0.63 10.5 + 0.89 14.3 + 0.86 
400 7.4+0.57 8.0+0.74 11.2+0.70 
500 6.84+0.58 6.15+0.78 7.9+0.67 

600 6.36+ 0,61 5.06+ 0.79 6.05+ 0.70 




Slcpci sid.dsv. 
o Pure CH 4 1-47 t 0.C3 
• 100% Af 1>25 i 0-09 
O 100% Ne l-l t 0-1 
A 100% He 0-2 Bi 0.03 


Log T^'K ^ 

Dependence 0 f ( T]/?)pyreCH 4 
on temperature. 


X 
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can be fitted by a power law, as in the earlier syatenie of 
the form (T^/^ ^CH -X^ where n=0.28 + 0.08 for CH^-He, 

1.1 + 0.1 for CH^-Xe and 1.25 + 0.09 for CH^-lr. The cr^ss 
sections ^ eff for spin-rotation interactions along ¥/ith the 
kinetic cross section sre tabulated in Table 11. 

Lalita and Bloom '' reported a value of 0.8 for the 
index n in CH^-He. They also have reported that T^/^ is not 
linear with composition for small concentrations of He and. 
assumed that these small departures do not appreciably affect 
the values of T^/p extrapolated to lOOfS He. Their values of 
T-j/^ were consistently smaller than the present values. It 
is possible that their results in the mixture were also 
influenced by the surface effects, as in pure CH^* IPig. 22 
shows that the present values of T^/^ in CH^-He mixtures fit 
very well the data reported by halita at room temperature 
for two different concentrations of He, where surface effects 
were believed to be eliminated. In the present results no non 
linear dependence of T-]/^ on the concentration of inert gas 
atom is found. Hence the extrapolated values of io 

100% inert gas atom can be taken as the real contribution to 

j, 

(T^/^ ) due to CH^-X collisions only where X could be He, le 


or Ar. 
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7.6 In to r pro ta. t lo n : 

As hef:.rc, defining A(J'M';JM) as the transition rate 
for the transitions JM J‘M‘ (CH^ moleqilQ) , it is given by 




A(J’M’;JI4) = 


2(J+1) 


V 


KK' 




-Cyp! 


j V^( t) I J‘K»M*) (J'K’M’i Y(t+T )) JIM)^dT 

( 5 . 10 ) 

The anisotropic potential Y_^(r, ©') is given by eq . 3.48. After 
carrying out the summations over K and K’ one obtains 


A(J'M';Jla) = [C(J3J'; (5.11) 


d(^'^tt«) = ^ 'i dr e X 

J 

— OO _____ 

(^b(r(t) Y*ye'(t)) b(r(t+t ) Y^^j^(ev(t+Z ))) 

( 5 . 12 ) 

with UOOjjt) —3(0). The function 3(0) is given by eq. 4.9. 

In case of CH^- Inert gas system, the index p which is determined 
by the angular variation of the anisotropic Interaction is given 
by 3. The radial dependence of the attractive part of the 
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anisotropic intortnolecular potential which is given by the 
index n = 7; the repulsive part of the anisotropic potential 
is assumed to vary as r"^^. The values of integrals '1(3,12), 
1(3,7) and 1(3,7,12) required for Ihe CH^-Inert gas atom 
potential are listed in Table 12 as a function of pe . 

How using eqs. 3.20, 3.21, 4.4 and 4.6, and the properties 
of Clebsch-Crordon coefficients, the life time of a molecule in 
a J state is obtained as. 



T92 -rc j(0) 
5145 - 


(5.13) 


and the relaxation time per unit density can be written as, 


f ^ 


5x10-” I 


1(3) 


(5.14) 


whe re 

1(3) = 1(3,12) + 1(3,7) -2 1(3,7,12) (5.15) 

7,7 Anisotropic Potential Parameters for CH^- Inert Gas Systems 
and Hyperpolarizability of Methane s 


The spin-rotation coupling constants were determined hy 
44 

Yt et, al for methane molecule from molecular beam experiments. 
The values are Cg_ = 10.4 KHz, G|j=18.5 KHz, The value of 0^^^ = 
(0^ + ( 4 / 45 ) 0?)"*^ = 138,6 MHz (Refer Table 3). The parameters 

3. CL 
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TA3I.E 1 2 

YIX^.:S3 OF I(T:, nO, I (p,n)< MI) I(p.n,nO M DILUHIE CLASSlCiil. 
GAS MIXIUIE FOR LMKIBB-JQNSS IS0TH3PIG POIEHTIALS 




1(3,7) 

1(3,7,12) 

0.06570 

0.2403 

0.0888 

0.1390 

0.07884 

0.1830 

0.0819 

0.1160 

0.10500 

0.1280 

0.0670 

0.0882 

0. 13140 

0.0969 

0.0612 

0.0730 

0.15750 

0.0780 

0.0560 

0,0620 

0.21000 

0.0568 

0.0493 

0.0495 

0,25640 

0.0460 

0.0458 

0.0435 

0.32050 

0.0380 

0.0433 

0.0380 

0.42733 

0.0297 

0.0405 

0.0325 
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and were obtained by fitting tbe expression 4.11 to the 

experiaental data of )„jj over the entire te'Bperature 

"*4 

region using least squares criterion. The hyperpolarizability A 
is given by (Refer Eq. 3.28), 


A = 


2 ^CH^-X ^2 


32 e 


( 5 . 16 ) 


CH^-X 


where a g is the polarizability of methane. The values a.j , Sg 
and A are tabulated in Table 13 along with the value obtained 
for CH^-Ar mixtures from infra-red spectral line broadening by 
Gray.^'^ 


V.S Discussion of the Results ; 

The value of the hyperpolarizability of me thane , A, 
obtained from CH^-lJe and CH^-Ar data agree well with the theore- 
tical value of 0..97 x 10~^^ cm'^ obtained from A = 4( 
a relation derived from a molecular model in which the tetrahedron 
is comprised of four atoms at a distance R^ from the center aid 
having polarizabilities a and a, along and at right angles 
to their bonds to the central atom^. The value obtained from 
CH^-He data is somewhat higher. Helium being the lighter atom 
than the other two, might be probing more of the repulsive part 
of the anisotropic potential than He and Ar at these temperatures. 
Hence it is likely that in this case the value of A is more 
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sensitive to the form of the repulsive potential. Further 

experiments at low temperature should clarify this point. 

The value of 1 obtained from CH,-Ar data may be compared 

57 

with the value obtained by Gray from the mean squared torque 
/ 2 \ 

m / on a CH^ molecule in a bath of Ar atoms as determined 

from Infra-red band moments. He has obtained a value of 
-32 4 

2.0 X 10 esu cm for A as compared to the present value of 
-32 4 

0.89 X 10 esu cm; Gray has the data only at 29 5° K and he 
could estimate only one of the two parameters involved, assuming 
the other to be zero. Hence it is possible that his value of 
A is overestimated. 


* * * 



CHAPOSR ¥I 


GONCLPSIONS 


It has been realized for quite sometime that the 
nuclear spin-lattice relaxation measurements can be inter- 
preted to obtain information on the intermolecular potentials, 
particularly the anisotropic part. measurements were 
indeed interpreted earlier In simple systems like hydrog’en, 
where the number of rotational states that are appireciably 
populated at temperatures of interest are restricted to one 
or two, to obtain precise information on the intermolecular 
potentials. So far such an analysis has not been done in poly- 
atomic molecules and their mixtures where the number of rotational 
levels are quite large. The present analysis is an attempt 
to obtain information on the intermolecular potentials, especially 
the anisotropic intermolecular potentials in polyatomic molecules 
and their mixtures from the experimental data of nuclear spin- 
lattice relaxation times, ‘T^s. 

Measurements of proton spin-lattice relaxation times 
were made in pure CH^, mixtures of CH^-C 02 and GH^- 

inert gas atom as a function of density and composition in 
the temperature region 300-600°K, using a phase coherent spin- 
echo spectrometer constructed for the purpose. The data 
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available in pure and v, 

. _ ^ ^4 2lso been interpreted 

0 oain similar information as in otb~T* « + 

in otbur systems studied. 

In the absence of a detail Ad 

detailed molecular theory for poly- 
atomic molecules to relata ^ J 

late the correlation time to the inter 
molecular potentials, it i. . T 

time Of the • " '""" " correlation 

Of ^ can be evaluated in terms 

theory ana Vea. oolllsion approrl.atlc„ 1 

Cl andS-p experimental data on pure gases CH., 

4 ®*^dSiP., two models nf ior^^- • . ^ 

have boon intermolecular potentials 

nave been considpr'A-^ 

i) » , - ' » ■ - appropriate anisotropic potentials, 

1 ) A hard sphere potential • ii'i 4 t 

case nt p ^ ’ ) A lennard-Jones potential. In 

case of hard sphere isotropic potential t>.P ■ ^ 

potential, the anisotropic part is 
assumed, to be desoT^ihAo^i -u - b .s 

described by an attractive term only and no 

repulsive part is considered m-hn 

. * necessary because when 

the rootrople potential is aeseri.ei .y a .ari sphere pciantial, 
oer am rntegrals in^l.ing anisotropic Intertoleoular functions 

become Inde.Penaent of temperature. Hence it is not r -a, * 

ficuue IX IS not possible to 

evaluate both the potential parameter^ repulsive as well aa 
attractive, of the anisotropic potential from the temperature 
aepenaenoc of T^/p alone. However, In lennara-Jones isotropic 
moael, in order to explain the observea temperature dependenoe 
of , it IS necessary to have a repulsive term also in the 
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anisotropic part alone witl an appropriate attractive term, 

' The scalar values of the octo|X 3 le moments of CH^ and 
obtained from the analysis of a hard sphere isotropic 
model are 20- 305fa less coaipared to the values from a lennard- 
Jones isotixjpic model. The oc to pole moments obtained from 
this model are in better agreement ?/ith the values reported 
from other techniques such as collision induced infrared 
absorption. In case of SiP^ the value of the octopole moment 
obtained from a hard sphere model is about 25^ higher compared 
to the value from a Lennard-Jones model. So far no value of 
Octopole moment of SiP^ is knovm from other techniques. 

The teroperature dependence of the proton spin relaxation 
data in CH^-COg and CH^- inert gas atom mixtures is 

different from a 3/2 lay; which is found to he the case in pure 
gases reported so far where spin-rotation interaction is the 
dominant mechanism of relaxation. These results have been 
interpreted as the failure of the hard sphere potential to 
adequately describe the isotropic potential. Hence a Lennard- 
Jones potential which is more realistic than a hard-sphere 
potential, is chosen to represent the isotropic part, with 
appropriate anisotropic potentials. The analysis yields ani- 
sotropic potential parameters and molecular quantities such 
as the octopole moment and hyperpolariz ability of CH^. The 
octopole moment of CH^ obtained from the analysis of 
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data is found to be 2.12 x 10 ssu cm® v/hen tiie qxiad.rupole moment 

P . 

of Ig is taken to be 1.5 x 10~ esu cm.This generally accented 
value of cuadrupole moinent of Ngr when corrected for hyperpolari- 
2 ability effects, is decreased by 20^, If the corrected value 
is used, the present anal 3 ''sis gives a value of 2.64 x 10”'^^®®'^ 
for the octopole moment of CH^. This value of 2.64 x 
is in good agreement with the value obtained from our analysis 
on pure CH^ where Lennard-Jones isotropic potential was used 
with an appropriate anisotropic potential; The octopole moment 
of CH^ obtained from the analysis of CH^-C 02 data is however 
smaller by a factor of two compared to the value from 
data. In case of CH^- inert gas atom the value 0.89 x 10~^^ cra^ 
that was obtained for the hyperpolarizabilitj?' A from CH^-So 
and CH^-Ar measurements is in good agreement with the theoretical 
estimate of 0.97 x 10 cm . However hyperpolarizahility 
obtained from CH^-He data is larger by a factor of t?’o. In case 
of CH^-C02 and CH^-He systems it is likely that the assumptions 
made regarding the radial term in the anisotropic potentials 
are not justified. 

It can be concluded that the present analysis of the 
experimental data on different systems studied indicate that 
the identification of the correlation time of the spin- 
rotation interaction with the average lifetime of a molecule 
in that J state is indeed a reasonable approximation. The 
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evaluation of the average lifetime under weak collision 
aprroxioation for appropriate anisotropic potentials f/ithin 
the frame work of liloom-Oppenheim theory, serves as a preli- 
rninar,'- estimate of the magnitude of the correlation time for 
spin-rotation interaction. It is also interesting to note 
that though the ratio of kinetic cross section to the effective 
cross section for spin-rotation interaction in all the systems 
considered indicate that the collisions are not really ’weak’ 
(mean number of collisions ranging from two to ten), ’the weak 
collision approximation’ can be successfully applied. 


* * * 
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NIEdERIGiL E7A1UA1I01 01 MlEGSil S 


General Remarks ; 

This Appendix presents a Rortran 17 program written 
for the numerical evaluation of integrals of type, 


I(p,n,n' ) 



It is noted that the integrals I(p,n) are evaluated from I 

j 

I(p,n,n') with n=n ' . The program includes two subroutines? 

i) A Gauss quadrature ii) A Bessel function routine. In the ; 

illustrative program presented here, + 32 point densities hav 

been used. The Bessel function routine contains 2 parts. Ro; 

smaller values of the arguement of J^^_i(xy) (say xy< 2) the i 

P+2 I 

following series has been used: | 


V 


(x) = (f) 




Mi 


m 


m=0 


m 


I [*" ( ) 


{|)2“ 


where ‘^ = p+-|-, Ror arguments greater than 2, the fractional 
order Bessel functions (order p) written in terms of elementa: 
Sin and Cosine functions have been used. The Bessel function; 
routines were tested by comparison with standard tables for 
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many selected values of arguments in all the orders (p) and 
were found to be accurate. Por the values of p,n,n* used in 
work, convergence in the x integration is generally achieved 
with an upper limit of x not exceeding 5, The corresponding 
limit in the y integration is around 30-50. The ’execution' 
presented in the example program displays the convergence 
achieved for 1(5,6,12) at two different temperatures. 
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IJOB PEG 122, TIME 

PAGES NAME 

S. RAJ AN 


DATE 


SIB JOB MAP 



SIBFTC MAIN 



PEG! 22 

FORTRAN 

SOURCE LIST 

ISN SOURCE 

STATEMENT 



0 

H 

O 

MAIN 


c 

SRI RAMAJAYAM 

1 


DIMENSION POINT (32), WEIGHT(32) 

2 


COI#iON / RAJ AN / POINT, WEI GET ' 

3 


READ 700,(P0INT(I),I=1,32) 

10 


READ 700,(WEIGET(I),I=1,32) 

15 

700 

F0RMAT(7F10.5) 

l6 

5 

READ 800,IP,BE,1N,BN,AL 

20 

800 

EORMATCt 5,4F10.5) ' 

21 


PRINT 900, IP, BE, AN, BN, AL 

22 

900 

FORMATC 4X, *!>=* ,1 3, 4X, ^^BE=* , FI 0. 5, 4X, *N=* , FI 0 . 5, 4X, * 

*,F10,5,*4X,*L=*,F10.5) 

23 


CALL FLUN( ) 

2if 


CALL GAUSS(BE,AN,BN,AL,IP) 

25 


GO TO 5 

26 

10 

STOP 

27 


EtID 


PHG122 IBFiAP ASSEIffiLY MAIN 

NO MESSAGES FOE ABOVE ASSEl^LY 

PHG122 FOETRAN SOUBCE LIST 

ISN SOURCE STATEflERT 


0 

illBFTC 

GAUSS 

1 


SUBROUTINE GAUSS(BE,AN,BN, AL,IP) 

2 


DIMENSION P0INT(32) ,WEIGHT(32) 

3 


GOMi.ON / RAJ AN / POINT,¥JEIGHT 

4 


DIMEN SION SUMT( 1 ) , SUMY ( 64) , SUMZ( 64) 

5 


XINT=3.0 

6 


SUMT(1)=0.0 

7 


YL=0.0 

TO 

23 

YL=YL+XINT 

1 1 


DO 35 1=1 ,64 

12 

35 

SUMY(I)=0.0 

14 


DO 43 ,64 

15 

45 

SUMZ(I)=0.0 

17 


XL=0.0 

20 


IXL=0 

21 

30 

XL=XL+XINT 

22 


IYY=0 

23 


IZZ=0 
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aif sssuM=: 0 .o 

25 DO 20 IY=1,32 

26 YP=POIHT(IY) 

27 5 YP1=YL+(YP-1.0)*]OTLV2.0 

30 DI=3.0*AL**2A.O/YP1**2 

31 D2=EXP(-D1) 

32 SUM! =0.0 

33 SUM2=0.0 

34 DO 15 IX=1,32 

35 XP=POINT(IX) 

36 10 XP1=XL+(XP-1.0)*XINT/2.0 

37 Xyi=XPl*YP1 

40 XGPT=XPl*»(-6) 

41 XGPT1=XGPT-1 .0 

42 GX=2.0*BE*XGPT1*XGPT 

43 IF(GX.LT.60.0) 6X=EXP(-GX) 

46 IF(GX.GE.60.0) GX=0.0 

51 T1=GX*BESSEL(r/1 ,IP)/XP1**(i\I’1-l . 5) 

52 T2=GX*BESSEL(XY1 ,IP)/XP1*^(BN-1 .5) 

53 SUM1=SUH1+?/EIGHT(IX)-='T1 

54 SUM2=SUM2+WEIGHT(IX)*T2 

55 XP=-XP 

56 IF(XP'.LT.O.O) GO TO 10 

61 15 COKTINIIE 

63 IYY=IYY+1 

64 IZZ=IZZ+1 

65 SUFiY(iyY) = SUIfY(lYY)+XINT/2.0*SUM1 

66 ,30LZ(IZZ)=SUMZ(IZZ)+XIHT/2.0*SUM2 

67 SSSUM=SSBUM+WEIGHT(IY)-^D2*SIJI4Y(IYY)*SWIZ(IZZ) 

70 YP=-YP 

71 IF(YP.LT.O.O) I® TO 5 

74 20 GOHTINUE 

76 TEPM=SSSm'4*XINT/2.0 

77 IaL=IXL+1 

100 SUKT(1XL)=SUMT(IXL)+TERK 

101 ppiifT 900, xl,yl,ixl,suf:t(ixl) 

102 900 F0Ri4AT(4X,E15.8,4X,E15.S,4X,I5,4X,E15.8) 


PHG1 22 

ISN SO UPC E 
103 
106 
107 

PHG122 


PORTRAfI SO BPCE LIST GAUSS 

STATEMENT 

IF(YL,LT.40.0) GO TO 25 

RETURN 

END 

FORTRAN PR0GRAJ1 GAUSS 


USER FUNCTION SUBPROGRAIU REFERENCES 


BESSEL 


PHG122 


IBMAP ASSEMBLY GAUSS 


NO MESSAGES FOR ABOVE ASSEJffiLY 



I^HG122 FORTRM SOURCE LIST 

ISN SOURCE STATEMHIT 
0 $IBFTC BESSEL 


1 


FUNCTION BESSEL(X,P) 

C 


GENERATE BESSEL FUNCTION OF ORDER P+1/2 

2 


INTEGER P 

3 


PI=3. 141 59625 

k 


IF(X.GT.2.0) GO TO 115 

7 


GAM=SQET(PI) 

10 


NP=P+1 

1 1 


DO 5 1=1, NP 

12 

5 

GAM= GAM* ( FLO Ar( I ) -0 . 5) 

14 


PNU=FL0AT(P)+0.5 

15 


SUM=1 .0/(2.0**PNU)/G/iM 

16 


P0WER=2.0**PNU 

17 


FACT=1 .0 

20 


PR0C=1 .0 

21 


XM2=1.0 

21 


M=1 

23 

105 

I=M 

24 


IF(M.GT.l) GO TO 6 

27 


P0WER=P0WEE*4.0 

30 


PRaD.=-PR0D 

31 


FAC T= FACT* FLOAT (M) 

32 


GAM=GAM* (PNU+FLO AT( M) ) 

33 


XI«I2=XM2*X*X 

34 


TERK=PR0 D* XM2/P0 WER/FAC T/GAM 

35 


GO TO 7 

36 

6 

Z=X*X/FL0AT(M)/(PNU+FL0AT(M))/4.0 

37 


TEEM=-TEEM*Z 

40 

7 

SUM=SUM+TEEM 

41 


IF(ABS(TEEM/SUM) .LT.1 .OE-10) GO TO 1 5 

44 


M=Mf 1 

45 


GO TO 105 

46 

15 

BESS}IL=SUM* ( X**PNU) 

47 


RETURN 

50 

115 

A?-= 945 0 * SIN ( X) /X* * 5 

51 


A2=945.0*C0S(X)/X**4 

52 


A3= 420 . 0* SIN ( X) /X* *3 

53 


A4=105.0*C0S(X)/X**2 

54 


A5=15.0*SIN(X)/X 

55 


A6=C0S(X) 

56 


BESSEL=A1-A2-A3+A4+A5-A6 

57 


BESSEL=SQRT(2.0/(PI*X))*BESSEL 

60 


RETURN 

61 


END 

PHG122 


IBMAP ASSEMBLY BESSEL 


NO MESSAGES FOR ^IBOVE ASSEMBLY 

IBLDR— JOB 000000 
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MEMO R|T 

SYSTEM, INCLUDING IOCS 
FILE BLOCK ORIGIN 

NUMBER OF FILES - 2 

1. S. FBIN 

2. S. FBOU 


MAP 

00000 THRU 12251 
12260 

12260 

12303 


OBJECT PROGRAM 12326 THRU 21641 


1. DECK' 

MAIN 

t * 

1 2326 

2. DECK' 

GAUSS 


12616 

3. DECK' 

BESSEL 

1 

13533 

4. SUBR' 

INSYFB 

f 

14274 

5. SUBR' 

OUSYFB 

1 

14333 

6. SUBR' 

POSTZ 


14364 

7. SUBR' 

CNSTNT 

T 

14736 

8. SUBR* 

FPR 

t 

14745 

9. SUBR' 

FRD 


14746 

10. SUBR' 

10 S 

t 

14747 

11. SUBR' 

RWD 

t 

1 5226 

12. SUBR' 

ECV 

r 

16402 

13. SUBR' 

FCV 

t 

16650 

14. SUBR' 

HCV 

t 

16742 

15. SUBR' 

ICV 

t 

17045 

16. SUBR' 

XCV 

t 

17063 

17. SUBR' 

INTJ 

t 

17103 

18. SUBR' 

FFC 

t 

1741? 

19. SUBR' 

FPT 

1 

20041 

20. SUBR' 

FLUN 

I 

20455 

21 . SUBR' 

XEM 

1 

20472 

22. SUBR' 

XP3 

f 

21037 

23. SUBR' 

XPN 


21 1 10 

24. SUBR' 

LOG 

t 

21216 

25. SUBR' 

SCN 

f 

21363 

26. SUBR' 

SQR 

t 

21542 


(^INSERTIONS OR DELETIONS MADE IN THIS DECK) 

INPUT-OUTPUT BUFFERS 
UBUSED CORE 

**«OBJECT PROGRAM IS BEING ENTERED INTO STORllGE AT 


76637 THRUE 
21642 THRU 
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P= 5 BE= 

0.30000000E 01 
0.30000000E 01 

O.3OOOOOOOE 01 

0.30000000E 01 
0.30000000E 01 
0.30000000E 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
0.30000000E 01 
0.30000000E 01 
0.30000000E 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 

p= 5 be= 
O.3OOOOOOOE 01 
O.30OOOOOOS 01 
O.3OOOOOOOE 01 
O.BOOOOOOOE 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 

O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
0 . 30000000 E 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 
O.3OOOOOOOE 01 


EXECUTION 


29200 

N= 

6.00000 


0.30000000E 

01 

1 

0.60000000E 

01 

1 

0.90000000E 

01 

1 

0 . 12000000 E 

01 

1 

O.I5OOOOOOE 

02 

1 

O.I8OOOOOOE 

02 

1 

0 . 21000000 E 

02 

1 

0 . 24000000 E 

02 

1 

O.27OOOOOOE 

02 

1 

O.3OOOOOOOE 

02 

1 

O.35OOOOOOE 

02 

1 

O.36OOOOOOE 

02 

1 

t 

O.39OOOOOOE 

02 

I 

O.42OOOOOOE 

02 

1 

35040 

N= 

6.000000 

O.3OOOOOOOE 

01 

1 

O.6OOOOOOOE 

01 

1 

O.9OOOOOOOE 

01 

1 

0 . 12000000 E 

02 

1 

O.I5OOOOOOE 

02 

1 

O.I8OOOOOOE 

02 

1 

0 . 21000000 E 

02 

1 

O.24OOOOOOE 

02 

1 

O.27OOOOOOE 

02 

1 

O.3OOOOOOOE 

02 

1 

O.35OOOOOOE 

02 

1 

0 . 36 '^ 00000 E 

02 

1 

O.39OOOOOOE 

02 

1 

0 . 42000000 E 

02 

1 


N 1 = 12. 00000 L= 0 .( 

0 . 751371 9 E- 0 if 

0.10739335-01 

0 . 25978 577 E -01 
0 . 28388 557E-01 
0.29102499E-01 
0 . 295201 55 E -01 
0.297751 55E/01 
0.29903062E-01 
0 . 299546 53E-01 
0. 29973829 E-01 
0.2998371 5E-01 
0 . 29991 106 E -01 
0.29996229E-01 

0 . 29998906 E -01 

NI= 12.00000 L = 
0 . 736 1031 3 E -04 
0 . 10247665 E -01 
0 . 23716505 E -01 
0 . 26060639 E -01 

0.26726083E-01 

0 . 27075 f.u 87 01 

0.27289895B-01 
0.27405759E-01 
0 . 27457 277 E -01 
0 . 27476596 E -01 

0.2748 4278 E-01 

O. 27488037 E-O 1 
0 . 27491 968 E -01 
0 . 27494564 E -01 



